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Abstract. We set up a model for the evolution of the 
galaxy luminosity function, taking advantage of recent 
work that brought in some better understanding of the 
mass function for gravitationally condensed objects. We 
add to this a simple model of star formation that repro- 
duces the behaviour of the Tully-Fisher relation, in order 
to attach a luminosity to a massive halo with a given veloc- 
ity dispersion. The physics of cooling of the gravitationally 
heated baryonic component allows us to distinguish halos 
that become groups or clusters from those that eventu- 
ally form galaxies (possibly within the former objects). 
With our new mass function and our new application of 
the cooling criteria - which motivated this paper - we get 
a satisfactory and natural cutoff at the bright end of the 
luminosity function, the needed flat slope for faint mag- 
nitudes and the correct trend in colors (brighter galaxies 
are redder) within the framework of the hierarchical clus- 
tering picture. This infirms earlier claims that the latter 
was inadequate to reproduce the former observations. We 
find the velocity dispersion to be a much better parameter 
than mass or radius to characterize galaxies. This model 
of the salient features that may describe galaxies allows 
one to discuss galaxy evolution as a function of redshift, 
in number as well as in luminosity. We find that bright 
galaxies form at z ~ 2 from mergers with a rather quiet 
evolution afterwards, whereas small galaxies are the result 
of a continuous merging process active up to the present 
epoch. The transition is found to occur at the observed 
transition between bright spirals and small dwarf ellipti- 
cals or irregulars. The galaxy luminosity was larger in a 
recent past for bright galaxies, as has been observed in 
the CFRS survey. This is because the mass of gas in a 
typical L* galaxy such as the Milky Way is a small frac- 
tion of the total baryonic mass and thus star formation 
is already slowing down. The evolution in number, which 
is quite well controlled in our model, agrees reasonably 
well with the counts as a function of apparent magnitude. 
The quasar multiplicity as a function of redshift is also 
discussed. 



Key words: Galaxies: luminosity function , evolution - 
quasars - large scale structure of Universe 

1. Introduction 

The formation and evolution of galaxies can be stud- 
ied through numerical simulations or analytic (or semi- 
analytic) models. In principle, the former could follow with 
great details the numerous processes (growth of density 
fluctuations, galaxy interactions, cooling processes, star 
formation histories, effects of supernovae and metallic- 
ity...) which govern galaxy formation. However numerical 
simulations are much less flexible than analytic means, 
and they also have to use physical parameterisations to 
model key processes like star formation, not only because 
the required dynamic range would exceed current compu- 
tational possibilities but also because these physical pro- 
cesses are not yet very well understood. Since the final re- 
sults depend strongly on such parameterisations (Navarro 
& White 1993), analytic or semi-analytic models provide 
an attractive alternative way to study galaxy formation, 
as they clearly show the relative importance of different 
processes and the influence of the available parameters, 
which would be more difficult to understand from a com- 
plex N-body and hydrodynamic simulation. 

The standard calculations follow the ideas of Press & 
Schechter (1974) (see also Schaeffer & Silk 1985), tak- 
ing advantage of the progress made originally by Bond 
et al.(1991), currently referred to as the "excursion set 
formalism" . The results obtained by such models (White 
& Frenk 1991; Kauffmann et al.1993; Cole et al.1994; and 
subsequent articles) show a reasonably good agreement 
with various observations which suggests that the main 
features of the standard scenario (galaxies arise from gas 
which is able to cool within dark matter halos formed by a 
hierarchical clustering process) are correct (or provide at 
least a good approximation). However, several problems 
are still unsolved: these studies cannot recover simultane- 
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ously the normalization of the Tully-Fisher relation and of 
the B-band luminosity function, the slope of the luminos- 
ity function at the faint end is too steep, massive and very 
bright galaxies usually arc not redder than fainter ones 
contrary to observations and the models often predict a 
cutoff at the bright end which is too smooth (there are 
too many very luminous objects). We present in this arti- 
cle an analytic model for galaxy formation and evolution, 
based on a specific description of density fluctuations in 
the highly non-linear regime and an original application of 
the cooling constraints, which provides a solution to these 
issues and allows detailed predictions for many physical 
galactic properties together with their rcdshift evolution. 

Thus, the motivations of this study are to: 

- take advantage of an improvement (Valageas & Scha- 
effer 1997, hereafter VS) over the usual Press-Schechter 
mass function (Press & Schcchter 1974), that originates 
in the understanding of the density fluctuations in the 
deeply non-linear regime (Schaeffer 1984; Balian & Scha- 
effer 1989; Bernardeau & Schaeffer 1991) and amounts to 
include the important subsequent non-linear evolution by 
counting directly the overdensities at the epoch of inter- 
est rather than in the primordial universe as implied by 
the Press-Schechter approximation. Since in this approach 
objects are defined by taking a snapshot of the actual den- 
sity field, in order to get the mass function of dark matter 
halos there is no need to follow the merging history of 
the various objects. This is done implicitely since the time 
evolution of the density field is already built-in within this 
non-linear hierarchical scaling model (and well verified by 
numerical simulations, see Bouchet et al.1991; Colombi et 
al.1997; Munshi et al.1998; Valageas et al. 1998a). Note 
however that the merging history of the dark matter ha- 
los enters, but only marginally, when the star formation 
history is considered. This is discussed in Appendix C. 

- implement in a natural way the cooling constraints. 
Indeed, the difference between galaxies and groups (that 
may have similar masses!) is that the former are able to 
cool rapidly (Rees & Ostrikcr 1977; Silk 1977; White & 
Rees 1978) whereas the latter are not. We shall argue that 
for a given velocity dispersion, galaxies are made from the 
largest possible baryonic patch that can cool within a few 
Hubble times at formation. By contrast, a larger patch 
which cannot cool rapidly in a uniform fashion (due to its 
lower overall density which translates into a larger cooling 
time) will become a group or a cluster of galaxies while 
several high density sub-units will cool and form distinct 
galaxies. Indeed, the non-linear scaling model predicts 
strong sub- clustering so that large clusters automatically 
contain numerous very high density spots that with our 
method we count as galaxies. For a galactic mass M our 
constraints translate into a well-defined mass-dependent 
condition on the galaxy radius R = R(M), or equivalently 
on its density contrast A = A(M), rather than having a 
universal (i.e. mass-independent) density contrast for all 



objects. Our approach is, in this respect, fundamentally 
different from the earlier ones. 

We apply this formulation to two descriptions of grav- 
itational clustering, based i) on a Press-Schechter-like ap- 
proach that uses linear theory and ii) on our non-linear 
hierarchical scaling model (VS). We show that in both 
cases it leads to a strong bright-end cutoff as required. As 
a result, we 

- show how the previous puzzles (normalizations of the 
Tully-Fisher relation and B-band luminosity function, flat 
slope for the latter, colors) can be solved within a hier- 
archical scenario. We can note that although relating the 
luminosity to the mass of the galaxy implies a model of 
star formation with a few adjustable parameters, observa- 
tional constraints (Tully-Fischer relation, gas/stars mass 
ratio,...) allow one to derive a mass-luminosity relation for 
galaxies independently of any consideration of the galaxy 
luminosity function. Hence the agreement with the ob- 
served luminosity function builds confidence in our model. 

- obtain detailed predictions for the redshift evolution 
of galaxies. Indeed, since there is nothing special about the 
present epoch the same considerations predict the evolu- 
tion of the galaxy mass and luminosity functions in the 
past (as well as in the future!). 

- take advantage of the possibilities of our analytic 
model to get scaling relations, analogous to the Tully- 
Fisher relation, between various physical quantities (cir- 
cular velocity - luminosity - mass - gas/star mass ratio - 
metallicity - star formation rate) over well-defined ranges 
of galaxy luminosity (for instance). Obtaining such laws, 
which requires the use of an analytic approach, provides 
a deep level of understanding and a clear presentation of 
the global trends implied by models based on the stan- 
dard hierarchical scenario, independently of the details of 
both the star formation process and gravitational cluster- 
ing. We think this transparancy and flexibility makes it 
worthwile to develop analytic models like ours (despite 
the simplifications they involve) in addition to numerical 
approaches. 

- present a model for galaxy formation within a global 
description of gravitational structures and astrophysical 
objects which provides a unified consistent model for very 
different phenomena: from clusters (Valageas & Schaeffer 
1998) to Lyman-a clouds (Valageas, Schaeffer & Silk 1998) 
and galaxies. 

The organization of the paper is as follows: in Sect. 2 
we discuss the determination of the mass function, with 
two models, one, that we call PS approach, based on an 
educated guess of which fluctuations identified in the early 
universe are going to eventually become galaxies, as pre- 
scribed by Press & Schechter (1974) but with a mass- 
dependent density contrast, and a second one, based on 
our previous work, using the same considerations applied 
directly to the non-linear density field at the epoch under 
consideration, that we call non-linear hierarchical scaling 
approach. In Sect. 3 we discuss the cooling condition, as 



Valageas & Schaeffer: The Mass and Luminosity Functions of Galaxies and their Evolution 



3 



well as our star formation model. In Sect. 4 we construct 
the galaxy luminosity function at the present epoch using 
the PS model as well as the non-linear scaling approach for 
a critical and an open universe. The evolution effects are 
examined in detail in Sect. 5. Finally, in Sect. 6, these con- 
siderations about the mass function and its time evolution 
are applied to derive the quasar distribution as a function 
of redshift, with an improvement over the Efstathiou & 
Rees (1988) model due to our new mass function that re- 
places the original one derived using the Press-Schechter 
model. In the last section we summarize the results of 
this paper and discuss its similarities and differences with 
other papers on the same subject (White & Frenk 1991; 
Blanchard et al.1992; Kauffmann et al.1993,1998; Cole et 
al.1994). 

The details of our model are presented in Appendix. In 
particular, we describe our star formation model in App.B 
and App.C. We discuss the values of our parameters in 
App.F and we present scaling relations in App.G. We use 
H a = 60 km/s/Mpc throughout this paper. 

2. The galaxy mass function 

Mass condensations may be characterized by two param- 
eters, such as for instance velocity dispersion and radius, 
or cquivalently mass and radius. As will be explained in 
Sect. 3.1, we define galaxies as halos of mass M with a ra- 
dius R{M) given by a condition provided by virialization 
and cooling constraints. This can be written as a mean 
density contrast A(M) within the radius R(M), so that 
we have a one parameter family of halos. The relation 
A(M) will also depend on redshift, but since all halos we 
shall consider will have already virialized we have for all 
z and O , and for any mass M: A(M) > 177. 

2.1. Press-Schechter approach (PS) 

We shall first devise a "Press-Schechter approach" (PS), in 
the sense that we wish to recognize in the early universe, 
when overdensities still grow according to linear theory, 
the density fluctuations which will eventually become ha- 
los of mass M, density contrast A(M,z), at the redshift 
of interest z. When this density contrast is constant and 
given by the virialization condition (obtained from the 
spherical collapse model) this is simply the usual Press- 
Schechter approximation. 

Let us first consider the case = 1. According to 
the spherical model, an overdensity characterized by the 
density contrast S(t) given by the linear theory, S(t) cx 



(1 



will slowly separate from the general expan- 



sion, its radius will reach a maximum R m at time t m 
and collapse at time 2 t m when 5 = S c0 ~ 1.68. We 
shall assume that in fact the halo virializes at this col- 
lapse time at the radius R m /2, as implied by energy con- 
servation and virial equilibrium if kinetic energy is neg- 



ligible at the turn-around. Hence at time 2 t m the den- 
sity contrast of the halo is A c0 ~ 177. Then we assume 
that the density of this object does not evolve signifi- 
cantly so that (1 + A(i)) oc (1 + z)~ 3 . Thus, a given 
halo of mass M, at redshift z, with a density contrast 
A(M,z) > A c o, collapsed at the redshift z c such that 
(1 + A(M, z)) = (1 + A c0 ) [(1 + z c )/(l + z)f and the den- 
sity contrast attached to this halo by the linear theory in 
the present universe is 5 c (z c ) = S c q (1 + z c ). If we note 
v = <5o/cto(M), where (Jq{M) is the amplitude of the den- 
sity fluctuations extrapolated to z — by linear theory as 
usual, we get: 



S c0 fl + A(M,z)\ 1/3 



a {M) 



1 + A, 



c() 



(1+*) 



(1) 



function of mass and redshift. Then, as long as 
A(M, z) > 1, we have 



n = i 



3 A(M, z) 



1/3 



10 cr (M) 



(1+*) 



(2) 



Then we consider that each halo is characterized by its 
parameter v so that the mass fraction in objects between 
v and v + dv is: 



ll{v) ± = Jl e S/* dv 

V V 7T 



(3) 



where we assumed gaussian initial fluctuations and we cor- 
rected by a factor 2 as in the traditional Press-Schechter 
prescription. Then the comoving multiplicity function of 
dark matter halos is: 



r){M) 



dM 
~W 



r){v) 



dv 



— e 1 dv 
M 



(4) 



where po = ~p(z = 0) is the mean density of the present 
universe. When the density contrast A(M, z) is constant 
and equal to A c0 the formulation (4) is exactly the usual 
Press-Schechter multiplicity function. In the case Qo < 
1 , A = 0, we have S(t) oc D{z) where D{z) is the growing 
mode of the linear theory and the density contrast of a 
halo which virializes at redshift z c is A c (z c ). Hence we 
now get: 



O < 1 , A = : v 



5 c {z c ) 
MM) 



where z c is defined by: 



(1 + A(M,z)) = (1 + A c {z c )) 



l + z c 
1 + z 



(5) 



(6) 



The multiplicity function of dark matter halos is still given 
by (4). Note that we have the normalization condition: 



Jo 



(7) 
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2.2. Non-linear hierarchical scaling model 

We shall now devise a second method which deals directly 
with the non-linear regime. We define for each halo the 
parameter x by: 



x(M,z) = 



1 + A(M, z) 
l[R{M,z),z\ 



(8) 



where 



L 



d 3 n d 3 



r-2 



V 2 



6(^1,^2) with 



V = 



is the average of the two-body correlation function 
£2(^1,7*2) over a spherical cell of radius R. Then we write 
for the multiplicity function of these halos at a given red- 
shift z (see VS): 



t I{m)—— — — x H(x) — 
M M w x 



(9) 



while the mass fraction in halos of mass between M and 
M + dM is: 



,„,.dM 2 tt / \ dx 



(10) 



The function H{x) is a universal function that depends 
only on the initial spectrum of fluctuations and that has 
to be taken from numerical simulations although its qual- 
itative behaviour is well-known: H(x) oc x^~ 2 for small x 
with uj ~ 0.5 and H(x) cx x UJs ~ 1 exp(— x/x*) for large x 
with uo s <~ —3/2 and x* ~ 10 to 20. Moreover, it satisfies: 



i; 



x 2 H (x) — = 1 
x 



(11) 



It has been shown in VS that this function H{x) is quite 
close to a similar scaling function h(x) that is obtained 
from the counts in cells. Bounds to estimate the differ- 
ence between these two functions are given in VS while 
numerical checks are presented in Valageas et al. (1998a). 
Hereafter we use the function h(x) in place of H(x). The 
correlation function £ that measures the non-linear fluc- 
tuations at scale R can be modelled in a way that follows 
very accurately the numerical simulations (see VS for more 
details). 

2.3. Comparison of the PS model with the non-linear hi- 
erarchical scaling approach 

These two formulations look very different, especially for 
evolution effects. The power at which redshift enters the 
expressions (4) and (9) is different even in the exponential. 
It happens, however, that the parameter v of the linear 
theory can be expressed in terms of the parameter x of 
the non-linear scaling theory. We will shortly summarize 
this calculation for Q = 1 . A detailed account can be found 
in VS. Using an analytic expression for the evolution of 



the correlation function, one can write the non-linear cor- 
relation in the stable clustering regime as a function of 
the extrapolation of the linear correlation function to the 
epoch under consideration: 



£(R,z) 



10 



R 



3a 

[l + l(R,z 



a(R L ,z) 



(12) 



R 



where a is a parameter close to unity. Using this relation 
one can obtain for a power-law initial spectrum of index 
n: 



1 5+2 

— X 6 

a 



(13) 



which shows (VS) that the expression (4) deduced from 
linear theory using the Press-Schechter ideas can be 
brought into a form quite similar to the non-linear scaling 
one (9) with: 



h PS (x) 



2 5 



7r 6a 



■ n 5±n 
-X 6 



' cxp 



v (5+n)/3 /(2a 2 ) (14) 



This holds in the highly non-linear regime, that is for 
£ > A c and A > A c . The definition of galaxies we shall 
use, which considers that galaxies are objects which have 
already virialized (and have since undergone a significant 
cooling) ensures that we always have A > A c for all galax- 
ies at any time. The condition £ > A c implies that (14) 
will hold until the redshift Zf ~ 3 such that £ ~ A c . This 
similar scaling as a function of the parameter x, for z < Zf, 
explains at the same time the similarity of the numerical 
estimates done with the two mass functions, but also their 
difference since the scaling functions H (x) and hps(x) are 
definitely not the same. The coefficient in the exponential 
is ~ 0.05 to 0.1 in the non-linear scaling case whereas it 
is ~ 0.5 in the PS approach, which implies a considerably 
faster fall-off at large x. The powers of x, typically larger 
than unity, entering the exponent in (14) further increase 
this difference. Thus, as compared to the non-linear scaling 
mass function the PS prescription overpredicts the number 
of intermediate objects (y <~ 1 or x ~ 1) and gives fewer 
extreme halos (very small or very massive objects), see 
VS for details. For z < Zf the non- linear scaling approach 
will give more numerous bright galaxies than the PS pre- 
scription, but fewer ones at z > Zf. A detailed comparison 
with numerical simulations of the various mass functions 
one can define within the framework of the non-linear scal- 
ing model, and with the PS prescription, is presented in 
Valageas et al. (1998a). 

3. Galaxy evolution 

The formalism introduced in the previous section only 
deals with the formation of dark halos through the ac- 
tion of gravity. To infer the mass function of galaxies we 
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need to precise when such a halo will form a galaxy, that is 
we have to express explicitely the physical processes which 
lead to the formation of a galaxy. 

3.1. Galaxy formation: cooling constraints 

Following the picture of the spherical model, an object 
with a small overdensity at early times will gradually sep- 
arate from the general expansion, reach a maximum ra- 
dius R m at time t m and finally collapse at time 2 t. m . One 
usually assumes that the real halo will virialize at the ra- 
dius i? m /2 at time p t m , with p = 2. We consider that 
once the halo of dark matter is virialized, the gas can cool 
and fall into the potential well as it is no longer thermally 
supported. However, at time (p + q) t m on average we as- 
sume that another larger halo containing the previous one 
will eventually form and will merge the initial overdensity 
with neighbouring ones. Indeed in the usual hierarchical 
scenarios, which we consider in this article, larger scales 
turn non-linear later so that small objects get embedded 
within increasingly large mass condensations. If the gas 
contained in the small halo has not cooled sufficiently, it 
will be distributed over the newly formed overdensity. In 
other words, the gas has only the time interval q t m to 
cool and jail into the initial halo in which it was embedded 
in order to allow the latter to retain its individuality and 
possibly form stars to become a galaxy. Fig. 1 displays this 
sequence of events. Of course cooling does not stop at time 
(p + q) t m . It may go till the present epoch, and the infall 
of baryons as well as star formation may still be active at 
much later time. 



log(R) 

turn-around 




larger scale 
collapse 



r ■ log(a) 

/ 

t m Ptm (P+q)t m 

Fig. 1. Collapse of an overdensity. Here R is the radius of the 
external shell, which turns around at tm and virializes at p t m . 
Matter condensations must virialize and cool in order to form 
a galaxy before they get embedded, at a time (p + q) t m on 
average, within a larger condensation. 

Thus, we write this cooling constraint as t coo i < q t m 
where t coo i is the cooling time. Hence the amount of gas 
which can cool and fall to the center of the halo, and pos- 



sibly form a disk and stars, is the mass of gas which is 
enclosed by the cooling radius R coo i where t coo i = qt m . 
However, for high redshifts or small temperatures, this ra- 
dius can be larger than the virialization radius R V i r , where 
the averaged density contrast A is equal to A c (given by 
p). In this case, the actual extension of the halo where the 
gas comes from is R V i r , since the surroundings have not 
collapsed yet. Hence we define the radius of the halos we 
consider by the constraints: 

{tcooi < qt m curve Ca 
and (15) 
p > A c p curve C v 

where p is the density of the universe at the redshift 
we consider. The parameter q should be of order unity 
(we shall take q = 2). Thus, the halos we consid- 
ered in the definitions (4) and (9) of the mass func- 
tion will be characterized by an external radius equal to 
R =Min(i? coo ;, R V i r ). This defines their extension and the 
function A(M) =Max(A coo j, A vir ) introduced in Sect. 2. 
As explained in Appendix A we do not need to assume 
a smooth power-law density profile for dark matter ha- 
los. Indeed, objects defined by R coo i are not defined as 
being the central core of radius R coo i within a larger halo 
of radius R V i r assumed to have a power-law density pro- 
file. We directly consider mass condensations satisfying 
the criteria (15) without reference to the properties of a 
possible larger object. The advantage of this procedure 
is that clusters automatically contain several galaxies (as 
they should) instead of only one central huge galaxy, see 
also the discussion below in Sect. 3. 2. The cooling time is 
given by: 

13S! p? e m p kT 1 

tcooi = - o o = ~ i coo ;( mean ha l° density) (16) 

s 2 ilb MpA(i ) s 

where A(T) is the cooling function (in erg.cm 3 .s _1 ). We 
use the cooling function given by Sutherland & Dopita 
(1993) for a gas with primordial abundances. The expres- 
sion (16) for s = 1 is the evaluation for the mean halo 
baryonic density. Cooling may be much more efficient since 
the local density increases within the halo and the baryon 
density gets higher during the collapse. This means that 
the parameter s which enters (16) should be larger than 
unity. A crude estimate for s using a mean power-law 
density profile (with no allowance for the cooling to be 
increased due to baryon concentration that would lead to 
even larger a value: such an increase has been shown to 
exist in numerical simulations, but is naturally limited, 
see Teyssier et al.1998) shows it may indeed reach several 
units. 

It is convenient to express (15) as a relation virial tem- 
perature - density, so that the halos we consider can be 
completely described by their temperature. Hence we can 
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write the curve p\ (T) corresponding to the first condition 
of the system (15) as: 



Ca : p\ — 



3tt 



flo 2>p^m p kT 



n b 2sqpA(T) 



(17) 



The curve C v corresponding to the second condition of 
(15), which characterizes just-virialized objects, is given 
by: 



C v : Pv = (1 + A c (z)) po (1 + zf 



(18) 



where A c (z) is the density contrast at the time of virial- 
ization, at the redshift z wc consider, and po (1 + z) 3 is 
the mean density of the universe at that redshift. If O = 1 
this density contrast is a constant: A c ~ 177. The general 
behaviour of these curves is shown in Fig. 2. The curve C 
which describes the system (15) is simply given by: 



Pc (T) = M a x[p A (T) 7 p v (T)} 



(19) 



This relation temperature-density defines implicitely the 
relation A(M, z) we introduced in Sect. 2 to determine the 
mass function of galaxies. Note that there are only two pa- 
rameters: p and (sq). From the spherical model we choose 
p = 2 as a natural value, while the product (sq) will be 
given by the luminosity of the Milky Way (which implies 
a well defined mass of baryons), with the constraint that 
both s and q should be of order of a few units. 

3.2. Galaxies versus clusters and groups 

The curves C\ and C v are shown in Fig. 2. Objects below 
either one of these curves violate one at least of the above 
conditions. Objects above these curves satisfy the crite- 
ria, but it is easy to see that slightly larger size objects 
(and whence lower density objects) also satisfy the crite- 
ria. Hence we define galaxies as being, for a given tem- 
perature (velocity dispersion), the largest object (hence 
the one with the lowest density) that satisfies the crite- 
ria (15). Thus the locus of galaxy formation, curve C, is 
the higher of the curves Ca and C v in the p-T diagram. 
Objects lying below this curve will be groups of galaxies 
or galaxy clusters. These objects do not cool in a uniform 
fashion within a Hubble time at formation, but may have 
cooled by now (groups) or are still cooling (X-ray clusters). 
The consideration of the latter in some sense justifies our 
cooling requirement for galaxies: we clearly see clusters 
cooling at the present epoch but they do not constitute a 
single galaxy. However, cooling is known to be necessary 
for galaxies to form and we would argue that the above 
condition is the most natural one to avoid the baryonic 
component of the following merging event to mix com- 
pletely. Groups and clusters contain some of the galaxies 
we are considering in the present paper and they can be 
described in a consistent way by means of the same meth- 
ods developed in VS and used here to describe the galaxy 



distribution. Their multiplicity, counting the same objects 
that the ones we deal with here but grouped differently, 
will be calculated elsewhere (Valageas & Schaeffer 1998). 
In particular, the normalization conditions (7) and (11) 
show that if we integrate from M — up to M = oo we 
recover the total mass of the universe whatever the curve 
A(M). The latter only describes how one can divide the 
matter content of the universe and different choices sim- 
ply correspond to different classifications: one can count 
a group as a single large object or as the assembly of sev- 
eral distinct galaxies which are considered as individual 
entities. In this article we adopt this latter point of view 
since we are interested in galaxies themselves. In fact, the 
galaxy mass function does not extend down to M = 
because halos with a low virial temperature T < 10 4 K 
do not cool (due to inefficient cooling as well as heating 
by the UV background). These patches of matter (which 
fill most of the volume of the universe with voids) form 
the Lyman-a forest clouds which we describe in details in 
Valageas et al.(1998b). 

iog(p ) 




log(T/l K) 



Fig. 2. Temperature-density diagram given by the cooling con- 
straint. Galaxies are defined as objects which cool within a 
Hubble time at birth. Objects which had no time to cool (be- 
cause of their low density, at a given temperature, which trans- 
lates into a long cooling time) are either galaxy groups (which 
have cooled by now) or hot X-ray clusters. They contain galax- 
ies located on C. 



The curve C\ does not depend on redshift while the 
curve C v does. At fixed temperature T, p v (T,z) as de- 
fined by C v increases with z along with the average density 
~p(z) of the universe. Whence the locus of galaxy forma- 
tion varies with z and at early times it becomes identical 
to C v . 

This variation with redshift has important conse- 
quences. Galaxies, that at a given epoch settle on the curve 
C v because larger objects of lower densities are not viri- 
alized, at a later epoch will lie above the curve C v and 
whence according to our criteria will be embedded within 
larger objects. This provides for these galaxies a contin- 
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uous merging process with obviously an evolution of the 
mass as well as the number of such objects (and indeed 
an associated star formation). On the other hand, galaxies 
that settle on the curve Ca at a later time will be limited 
by the same condition on density: larger objects of lower 
density will not satisfy the above conditions and thus may 
be groups or galaxy clusters, but cannot become a larger 
galaxy into which the galaxy we consider gets embedded. 
So in the framework of our model these galaxies will no 
longer evolve neither in mass nor in number. We call these 
galaxies "isolated galaxies" . Thus we have a population of 
smaller galaxies that evolve by mergers and a population 
of larger galaxies that bear only internal evolution. At 
early times all objects are continuously merging since the 
only limiting curve at high z is C v . As time goes on, the 
larger objects get limited by cooling, settling on Ca- So 
all these now "quiet galaxies" have an early phase with 
strong merging processes. Once this early phase is over, 
the now quietly evolving galaxy may for instance form a 
disk. Within this picture, we do not take into account the 
strong merging processes that take place in the very dense 
cores of clusters. Such additional interactions may suffi- 
ciently disturb (Balland, Silk, Schaeffer 1998) the galaxies 
we model here to change them into ellipticals, a plausi- 
ble explanation of the density morphology relation. This 
will be accounted for in a subsequent paper about galaxy 
clusters (Valageas & Schaeffer 1998). 

Note that in our model all galaxies have roughly the 
same age (of the order of the age of the universe) from 
the star formation viewpoint, even if they did not exist as 
distinct objects along this whole period. Indeed, present 
galaxies are the result of the merging of many smaller and 
older sub-units where star formation was already active at 
high z. In fact, contrary to the usual statement, although 
we work within the framework of the standard hierarchi- 
cal clustering scenario massive galaxies look slightly older 
than small ones (see App.C). Because of their more effi- 
cient star formation process (due to their higher density 
and virial temperature, see App.B) massive and bright 
galaxies have a redder and older stellar population. 



3. 3. Star formation 

Our scenario describing the history of star formation is 
rather standard, and kept as simple as possible. Apart 
from a small component (10%) of dark baryonic matter 
(brown dwarfs, planets) that plays a negligible role in the 
stellar evolution history, we consider four components: 

- short lived stars, of total mass M s h, that will be re- 
cycled. 

- long lived stars, of total mass M; D , that will not be 
recycled. 

- a central gaseous component, of total mass M gc , at 
the sites of star formation, that is deplenished by star 
formation and ejection by supernova winds, replenished 



by infall from a diffuse gaseous component located in the 
dark halo potential well. 

- a diffuse gaseous component, of total mass M g h, 
deplenished by infall and replenished by the supernova 
winds. 

The star formation rate dM s /dt is proportional to the 
mass of central gas M gc with a time-scale t c proportional 
to the dynamical time-scale r^: 



dM s 
dt 



M, 



and t c cx Td 



(20) 



The mass of gas heated and ejected by supernovae out of 
the central parts of the galaxy is proportional to the star 
formation rate and decreases for deep potential wells: 



dM g 
dt 



SN 



To MgC 
T T c 



(21) 



where T is a constant and T the halo virial temperature, 
see (B2). Finally, the infall of gas from M g h to M gc occurs 
on a dynamical time-scale and not on the cooling time- 
scale since t coo i < Td because of (15): 



dMgA = Mgh 
dt J IN r d 

Since Td ~ 1/ \[Q~P ~ t m we shall use: 



Td = (3d t„ 



Tc — Pc tn 



(22) 



(23) 



where (3d and (3 C are parameters of order unity. We as- 
sumed t c oc Td, as it may describe gravitational instabili- 
ties within galaxies as well as the influence of neighbours. 
Moreover, the time-scale r c disappears for faint galaxies, 
as we shall see below, and our model would still be ap- 
proximately valid even for bright galaxies if one has in fact 
t c < T,;, since it would simply correspond to a change of 
[3d (then the system is governed by the longest time-scale 
among r c and Td). This means that other forms for r c , for 
instance t c cx l/p or r c cx l/p g (where p g is the gas den- 
sity in the core) would give similar results, as we checked 
numerically. 

The details of this model of star formation are pre- 
sented in App.B and App.C. We solve the corresponding 
equations (Bl) numerically but we also give there approx- 
imate analytic solutions that allow one to follow the be- 
haviour of these various components as a function of time 
and galactic mass. We can stress here that the fraction 
of non-luminous stars we adopt (10%) is consistent with 
observations (see discussion in Mora et al.1998) while the 
values used in some other studies (50% in Kauffmann et 
al.1993; 63% in Cole et al.1994 for instance) are much too 
large. Using a smaller fraction of non-luminous stellar- 
like objects however would increase the luminosity of the 
galaxies obtained in the latter models and would lead to 
stronger discrepancies with the observed luminosity func- 
tion. 
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We also describe the evolution of the metallicity (El) 
with three different values: Z s for stars, Z c for the central 
gaseous component and Zh for the diffuse gaseous compo- 
nent. 

We summarize here the free parameters of our model: 

- (sq): proportionality factor in the cooling constraint 
(17). 

- (3 C : proportionality factor for the star formation time- 
scale (23). 

- Pd- proportionality factor for the definition of the 
dynamical time-scale (23). 

- T : supernovae efficiency, (B3) and (21). 

- rjd = 10%: fraction of non-luminous objects (brown 
dwarfs) (D8). 

- "yield" y — 0.019: mass of metals produced per mass 
of stars through SNII, see (El). 

The parameters q, (3 C and (3d must be of order unity 
and they are constrained by the luminosity, mass, gas/star 
mass ratio of the Milky Way, while T is constrained by 
the supernovae rate of similar galaxies. The yield y is given 
by the metallicity of the solar neighbourhood. Thus all pa- 
rameters of our model are adjusted on other observations 
than those related to the luminosity function we seek to re- 
produce. Note moreover that for bright and massive galax- 
ies (T To) To is irrelevant: the effect of supernovae is 
negligible since the potential well is very deep so that gas 
cannot be ejected efficiently For faint galaxies (T <C To) 
(3 C is irrelevant because a quasi-stationary regime sets in 
very quickly where the galaxy evolution is governed by 
the equilibrium between the ejection of central gas by su- 
pernovae and the infall from the diffuse component (see 
App.B.2) so that the star formation rate is given by: 



dM s 
dt 



TO 



with t ~ 1 + 



Td 



(24) 



where M g is the total gas mass. For massive galaxies the 
derivation leading to this relation is no longer valid but 
(24) still gives a reasonable description of the star forma- 
tion rate because in our model Td ~ t c so that for T 3> Tq 
there is only one time-scale which is also given by To since 
in this case To — Td- The value of the main parameters 
(sq) and (3d is further discussed in App.F. 

4. Numerical results 

Now we can compute numerically the galaxy luminosity 
function, as denned in the previous section, as well as the 
physical characteristics of the halos we consider. 

4.1. fl = 1 

4.1.1. Halo properties 

In the case Q — 1, Fig. 3 shows the relation temperature- 
density which corresponds to the cooling and virialization 
constraints (Fig. 2), at redshift z = 0. We see that the 



typical mass and radius of the halos we get are close to the 
observed values for galaxies in the present universe. For 
T ~ 10 6 K the cooling constraint translates into a nearly 
constant mass M ~ 10 12 M Q , while at large masses, or 
large temperatures T > 10 7 K, it translates into a fixed 
radius R ~ 100 kpc. Indeed, the curve C\, see (17), can 
also be written: 



R = (sq) 



fit / ( 7 -i)MCQ 2 

fi V 2£ 2 M 4 m 3 kT 



(25) 



At high temperatures, we can write A(T) ~ Kq^/T/Tq 
with T = 10 7 K and A = 10~ 23 erg.cm 3 .s _1 . This gives 
a constant radius R: 



(26) 



The fact that we know the radius of the halos we consider 
to be at least 60 kpc (because a few rotation curves mea- 
sured in large spiral galaxies remain flat at least up to this 
radius) requires that (sq) > 7. However, as we explain in 
App.F the product (sq) is determined independently by 
the luminosity of the Milky Way. This leads to (sq) ~ 12 
in the case fi = 1, fit — 0.04. Hence both constraints can 
be satisfied simultaneously. Since we use p = 2 as given by 
the simple spherical collapse picture and we expect q ~ p 
and we noticed in Sect. 3.1 that s > 1, we choose q = 2 
and s — 6. This gives R — 102 kpc in the limit of very 
large virial temperatures. 

Galaxies with small virial temperature T < 10 6 K and 
rotation velocity V c < 100 km/s are irregular or (dwarf-) 
elliptic galaxies, dominated by their continuous merging 
history, while massive galaxies with V c > 100 km/s are 
disk galaxies which have evolved through a much calmer 
history since they reached the curve C\. Indeed, they re- 
main unchanged, or slowly accrete some mass if cluster- 
ing is not exactly stable, so that a disk can form. This 
boundary V c = 100 km/s (which corresponds to a B-band 
magnitude Mb = —18) is indeed close to the observed 
transition between spirals and (dwarf-) ellipticals or irreg- 
ulars (Sandage et al.1985). 

We can see that in the near future (z < 0) two halos 
of the same mass or the same radius can have different 
temperatures and densities. This implies that, within the 
framework of this model, neither the mass nor the radius 
are good variables to describe the halos we consider, while 
the velocity dispersion (or the temperature) is. Note that 
the latter is very close to our variable x that we think is 
the true parameter characteristic of mass condensates. 

We notice that light and faint galaxies, which cor- 
respond to small temperatures or velocity dispersions, 
are young objects characterized by a density contrast 
A ~ 200. Indeed, these small galaxies are on the curve 
C v in the cooling diagram (see Fig. 3). This may look sur- 
prising, since the observed density contrasts of these ob- 
jects are usually very large. However, this is simply ex- 
plained by the fact that the observed effective luminous 
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Fig. 3. Temperature-density diagram given by the cooling 
and virialization constraints in the case f2 = 1, z — 0. The 
short-dashed lines correspond to halos of fixed mass, on the 
left, and halos of fixed radius, on the right. We note p c o the 
critical density at z = 0. 

radius Ri um is much smaller than the radius we consider 
which describes the actual size of the dark matter halo 
from where the gas content of the galaxy originated. This 
large difference in radius translates into the density con- 
trast we obtain. The small size of the observed luminous 
radius of these objects, that is not modelled here, can be 
explained by several effects. For instance a large fraction 
of gas is ejected, since the potential well is too weak to 
retain baryons as efficiently as in large galaxies. Also star 
formation is not very efficient so that it cannot take place 
very far from the center of the galaxy because the density 
gets quickly too small. 

We do not attempt to model here the radius and the 
mass within the luminous part. This has already been done 
in a similar context (Bernardeau & Schaeffer 1991). It was 
seen there by phenomenologically modelling the baryon 
squeezing due to cooling, that the observed behaviour of 
the luminous radius calls for a dark matter halo radius 
that behaves nearly as M 1 / 3 for the smaller masses (i.e. 
a constant density contrast) but for an almost constant 
dark halo radius for the large masses (due to a change 
in the M/L ratio in these two regimes). These findings, 
indeed, encouraged us to undertake the present work. For 
instance, if Ri um cx L 5 and the density profile of galactic 
halos is a power-law p(r) cx r -7 (7 ~ 1.8), then for small 
galaxies on C v where L cx V 3 the density contrast within 
the luminous radius scales as (1 + A); um cx 

L-t/ 6 w hilc for 

bright galaxies on C\ which have a nearly constant dark 
matter radius we get (1 + A)i um cx £~( 3 t~ 4 )/ 6 . Hence the 
density within the luminous radius decreases for brighter 
galaxies even though the density within the much larger 
dark matter halo increases. However, a precise estimate 
of these relations would require a detailed model of the 
behaviour of the gas after cooling which we do not consider 
in this article. 



4.1.2. Gas/star mass ratio 

The variation of the gas/star mass ratio, which is closely 
related to the ratio (star formation time-scale) / (age of the 
galaxy), see (D4), is shown in Fig. 4 as a function of the 
galaxy circular velocity. Small V c galaxies which have a 
small star formation rate Tq 1 , because of their low density 
and virial temperature, have a high gas/star mass ratio 
M g /M s 3> 1. On the other hand, high temperatures at 
the external radius of the halo, or large V c , correspond 
to high densities. Hence these galaxies have undergone 
a very efficient star formation, since their star formation 
rate Tq 1 is large. As a consequence, they have already 
lost most of their gas content, which turned into stars, 
as we can see in Fig. 4, and M g /M s <C 1. As expected 
on general grounds, see (D5), we verify that M g /M s <~ 
1, corresponding to a galaxy similar to the Milky Way, 
implies r ~ t (we have M gc /M s = 0.37 and M g /M s ~ 
1 since there is some gas in the halo: M g h — 2 M gc ~ 
0.008 M). This value of the gas/star mass ratio leads us 
to choose for our last parameters fi c = 2, [3d = 3 and 
Tq = 3 10 6 K. We can see in Fig. 4 that while most of the 
gas is in the dense component M gc for galaxies brighter 
than the Milky Way (with a larger rotation velocity) , it is 
mostly in the diffuse phase M g h for faint galaxies, as we 
can see in (B5), because supernovae are very efficient in 
these weak gravitational potentials and eject the gas out 
of the star-forming regions. In particular, from (B5), (B9) 
and (BIO) we see that for faint galaxies M gc ~ M s while 
M gc ~ M g for bright galaxies. 




log(Vj [ km/s ] 

Fig. 4. Ratios of the total mass of gas M g (solid line) and of 
the mass of star-forming gas M gc (dashed line) over the mass 
M s of luminous stars as a function of the circular velocity V c . 

One may distinguish three different regimes (App.G) 
in which approximate power-law relations can be derived. 

1) Very faint galaxies, located on C v : A ~ 200 and 
T < 10 5 ' 5 K, with a constant density contrast. 

2) Faint galaxies, located on C\ ■ R — 120 kpc and 
lO 5 5 K < T < lO 6 5 K, with a nearly constant radius. 
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3) Bright galaxies located oiiCa : R — 100 kpc and 
T > 10 6 - 5 K, which have nearly exhausted all their gas. 

Using these analytic estimates (App.G) it is readily 
seen that the gas/star mass ratio is a steep function of 
the circular velocity (M g /M s oc V~ 2 to M g /M s oc V~ 3 
and finally it follows an exponential decline). We can see 
very distinctly in Fig. 4 these three regimes. 

It is important to realize that these relations should 
be approximatively valid for any model of star formation. 
The relation (D5) is always correct for t < t , where t is 
the star formation time-scale. Indeed, for t < tq we have 
dM s /dt ~ M 9 /to by definition of t , hence M s ~t/r M g . 
Since the Milky Way is characterized by M g /M s ~ 1, 
fainter galaxies (which have a higher gas/star mass ratio) 
verify t < t , M g ~ M g0 and M s ~ t/ro M g0 , while more 
luminous galaxies satisfy t > t and M s ~ M g0 . This 
means that the relation V c — L for these bright galax- 
ies, if the stellar mass/luminosity ratio M s /L is roughly 
constant, is directly given by the physical characteristics 
of the underlying halos - since L oc M s a M - hence 
by the cooling curve, and not by the details of star for- 
mation processes. Hence the slope of the corresponding 
Tully-Fisher relation (which is a bit shallow) could only 
be modified by a change of the definition of galactic halos 
(i.e. using another constraint than the cooling criterion) 
or by a gas/dark matter mass ratio which would vary with 
the galactic characteristics (one may argue that deep po- 
tential wells could gain some gas from surrounding small 
halos). On the contrary, the luminosity of faint galaxies 
L oc t/ro M depends strongly on the star formation rate. 
Using p = A c p and t~tn (where ta is the age of the uni- 
verse), we have L oc U 3 /to- Hence the dependence of To 
on V is constrained by the observed Tully-Fisher relation. 
In fact in our model the relation we obtain in this case 
is To oc V~ 2 but it leads to results which are still consis- 
tent with observations. Thus, the high uncertainty which 
lies in any star formation prescription (since this process 
is rather poorly known) is greatly reduced by these gen- 
eral properties and the observed Tully-Fisher relation. It 
is also apparent through these considerations (and is con- 
firmed by a simple calculation) that the use of a constant 
star formation time-scale (r = To for all times) would not 
change these results (at z = 0), except for small details. 
As a consequence, we can reasonably expect our results to 
be fairly general and robust. 

4.1.3. Luminosity 

Galaxies characterized by high temperatures and large V c 
also correspond to the most luminous galaxies, and their 
luminosity is mainly due to small stars which have a long 
life-time, since there is not much gas left to create new 
generations of massive short-lived stars, as it appears in 
Fig. 5. We can see in the lower panel that the Tully-Fisher 
relation is approximately satisfied, although the slope of 



the relation V c — L gets shallower at high luminosities. For 
B-band luminosities, observations give L oc V 2 - 7 (Kraan- 
Korteweg et al. 1988; Pierce & Tully 1988), while the in- 
frared relation is steeper: L oc U c 3 5 (Pierce & Tully 1992). 
In our model, these relations arc obtained (App.G) as a 
smooth transition from L oc V 5 at the very faint end to 
L oc V 2 at the bright end. 



o.u 




log(V ) [ km/s ] 



Fig. 5. Upper figure: ratio of the total mass of luminous stars 
M s over the galaxy luminosity L. Lower figure: luminosity L of 
the galaxy as a function of V c . The thin dotted line corresponds 
to the power-law L oc U c 3 normalized to the Milky Way. The 
short-dashed curve is the luminosity due to long-lived small 
stars (Li Mi ) while the dot-dashed curve is the luminosity 
due to short-lived massive stars (L s hM s h). 

Hence the relation V c — L we obtain has a strong 
slope close to the usual Tully-Fisher relation, and gets 
shallower for galaxies more luminous than the Milky Way. 
It is interesting to note that Persic & Salucci (1991) found 
a similar decrease of the slope of the Tully-Fisher relation 
for luminous galaxies in observations, although in the in- 
frared H-band. We can note in Fig. 5 that for faint galaxies 
the slope is close to 3 and not 5 as in the approximation 
obtained in App.G because the ratio M s /L is not exactly 
constant and increases for large circular velocities. 

As we can see from (Dl) and (D4), this variation of the 
mass/luminosity ratio M s /L is due to the fact that large 
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and bright galaxies (also characterized by a large tempera- 
ture or circular velocity) , which have a high star formation 
rate Tq 1 , have already consumed most of their initial gas 
content. Hence, their present star formation rate is rela- 
tively small (compared to their past history), and their 
stellar population presents an increasing proportion of old 
stars, with a long life-time, which were created during the 
whole life of the galaxy. Thus, as we can see in Fig. 5 the 
luminosity in the form of short-lived massive stars is much 
smaller than the contribution of long-lived small stars. On 
the contrary, small galaxies which still have a relatively 
high star formation rate (as compared to their past) be- 
cause of their large gas content show an important contri- 
bution from massive stars which are created at the current 
epoch. In fact, both classes of stars give roughly the same 
luminosity. This is due to the fact that the luminosity 
of the galaxy is dominated by stars of intermediate mass 
<~ 1 Mq, as described in App.D. This variation of 
the mass/luminosity ratio shows that very massive galax- 
ies should be redder than small ones, which is consistent 
with observations (Lilly et al.1991; Metcalfe et al.1991). 
Moreover, this would be enhanced by metallicity effects, 
since bright galaxies are also the most metal-rich, which 
we recover in our model as we shall see now. 

4.1.4. Metallicity 

Next, we can also consider the metallicities Zh (diffuse 
gas), Z c (dense gas) and Z s (stars), which we introduce in 
App.E. They are displayed in Fig. 6, as a function of the 
B-band luminosity of the galaxy. For faint galaxies we are 
in the "stationary" regime where Zh <C Z c and Z s ~ Z c . 
The variation of the metallicities Z c and Z s , which are 
those available to observations, agrees with the data. We 
chose the value of the yield y so that we get Z c = Zq for 
the Milky Way. This gives y = 0.019 ~ Z Q . 

As usual, we can consider the 3 regimes we intro- 
duced previously to determine approximate relations for 
the metallicity (App.G). For galaxies fainter than the 
Milky Way we have Z c cx Z s cx L 0A , or L cx Z„. Galaxies 
more luminous than the Milky Way have a constant metal- 
licity Z s ~y while Z c cx L 0,5 . Note however that the yield 
may vary with the characteristics of the galaxy, together 
with the IMF. Nevertheless, we see that our model agrees 
with observations over 10 magnitudes in Mb- 

4.1.5. Stellar history 

Finally, we can consider (rather crudely) the mean mor- 
phological properties of the galaxies we obtain in our 
model as a function of their B band magnitude Mb- We 
define an approximate disk/bulge luminosity ratio by: 

_ L sh M sh + L l0 (Mi (t H ) - M io (p t m )) 
Li Mio(pt m ) 

when tn > p t m (tn is the Hubble time at the considered 
redshift and p t m the virialization time of the galaxy). 




- IO - 15 -20 

Fig. 6. Metallicity Z as a function of the B-band magni- 
tude Mb- The solid line corresponds to the gas in the core 
(Z c ), the short-dashed line to stellar metallicities (Z s ) and the 
dot-dashed line to the gas in the halo (Zh)- The data points 
are observations from Zaritsky et al.(1994). 

Thus, we evaluate the disk luminosity as the contribution 
from stars formed after the galaxy reached the cooling 
curve C\, and the bulge luminosity as the contribution 
from earlier stars. For galaxies which have not reached 
Ca yet we use L^/b = since there is no disk. This ap- 
proximation is only very crude, as the disk may not form 
exactly at the virialization time p t m , and some stars prob- 
ably form in the bulge after this date. The variation of this 
ratio L d /b with Mb is shown in Fig. 7. 




O I ' 1 ' ' ' 1 ' ' ' 1 ' ' 1 

-22 -20 -18 - 16 
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Fig. 7. The variation of the disk/bulge luminosity ratio L D / B 
with the B band magnitude M B - 

Faint galaxies Mb > — 18 have not reached Ca yet, 
hence they have no disk and are elliptical or irregular 
galaxies. For Mb < —18 the disk/bulge luminous ra- 
tio first increases with luminosity, since brighter galax- 
ies (which also have a higher density, virial temperature 
and a larger mass) have been on Ca since a higher red- 
shift. However, this ratio declines for very luminous galax- 
ies Mb < —21, which correspond to large densities and 
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deep potential wells, hence to very efficient star forma- 
tion. As a consequence, these galaxies transformed a large 
part of their initial gas content into stars at very high 
redshift, while they were still divided into several sub- 
units, before they reached the curve Ca which marked 
the end of their merging phase of formation. Hence, the 
mass of stars formed since this latter epoch is increas- 
ingly small as compared to the stellar population formed 
during the merging phase, as the galaxy parameter x in- 
creases. Thus, it appears naturally in our model that spiral 
galaxies only correspond to an intermediate range of lu- 
minosities, —22 < Mb < — 18, while brighter and fainter 
galaxies should be ellipticals or irregulars. Hence we get 
old bright ellipticals as a straightforward outcome. This is 
quite similar to the observed dependence of the dominant 
Hubble type on the luminosity (Sandage et al.1985). How- 
ever, the interactions between galaxies in clusters would 
certainly add several effects which we did not take into ac- 
count explicitely but could have important consequences 
on the mean luminosity-morphological type relation. 



4.1.6. Milky Way 

The Milky Way corresponds to V c ~ 220 km/s, and 
in our model to R ~ 140 kpc, M ~ 2.6 10 12 M Q , 
A ~ 2200, T ~ 3.6 10 6 K , L ~ 2 1O 1O L , M B ^ -20.4, 
M gas ~ 5.5 10 10 M Q and M gc /M s ~ 0.37. Its present 
supernovae rate is Rsn — 0.025 year -1 which is con- 
sistent with observations: van den Bergh & Mc Clure 
(1994) find R SN = 0.021 - 0.024 year" 1 . The above 
value of Mg C /M s is consistent with the values a gas /<7 s of 
the disk surface densities measured in the solar neighbor- 



hood: a gas /(T s = 0.27 (Bienayme et al.1987), 0.47 (Kui- 
jken & Gilmore 1989), 0.5 (Gould et al.1996). It implies 
a gas/(gas+stars) mass ratio of <~ 0.3 which is closer to 
the actual observations than the value ~ 0.1 traditionally 
taken in many calculations of the chemical evolution in the 
solar neighborhood. The Milky Way also corresponds to 
a star formation rate dM s /dt ~ 5M Q /year, which agrees 
with usual estimates. Note that the present ratio (stellar 
mass)/(age of the universe) M s /t ~ 4 (in M©/year) im- 
plies that any model with the correct luminosity -hence 
the correct stellar mass- will give a present-day star for- 
mation rate of the right magnitude. 



where 

I = (flak 2 )- 1 , a = 1.7 I , p = 9 l 3/2 , 7 = I 2 

and A is a normalization constant such that as = 0.5. 
Fig. 8 shows the linear correlation function <j(R) 2 and the 
evolved non- linear correlation function £(i?) at redshift 
z = 0. The calculation of the non-linear correlation func- 
tion £ from its linear counterpart a 2 , using the spherical 
model normalized by the numerical calculation performed 
by Jain et al.(1995), is detailed in VS. 




log(R) 



Fig. 8. The evolved correlation function £(R, z) at z = 
(solid line), and the linear extrapolation u(R) 2 (dashed line), 
in the case f2 = 1 and P(k) is a CDM power-spectrum. The 
dot-dashed line is the power-law (_R/i? )~ 7 with 7 = 1.8 



We can see that the evolved non-linear correlation 
function is close to a power-law £(R) oc i?~ 7 , with 7 = 1.8, 
in the range 50 kpc to 1 Mpc, which covers all the values 
that the radius of the halos we consider can take. As we 
saw in Sect. 2, we can now get the luminosity function of 
galaxies, provided we know the function h(x) we intro- 
duced in that section. We choose the function h(x) given 
by Bouchct ct al.(1991): 



h(x) = a 



with: 



1-W 



T{u) {l+bx) c 



-x I X* 



(29) 



4.1.7. = 1, CDM 

Now, to get the mass function or luminosity function of 
galaxies, we need the value of the initial power-spectrum 
-P(fc), or the correlation functions cr(R) 2 and £,(R). We 
first consider the case where P(k) is a CDM-like power- 
spectrum. More precisely, following Davis et al.(1985) we 
use: 

P(k) = A k (1 + ak + (3k 1 - 5 + -(k 2 y 2 Mpc 3 (28) 



a =1.8 , 6 = 3.6 , c = 0.8 , w = 0.4 , = 12.5 

The galaxy luminosity function we get in this way is 
shown in Fig. 9 as a function of the B-band magnitude 
(M B = 5.48 - 2.51og(L/L Q )). The short-dashed curve is 
the prediction of the PS approach. 

We explained in Sect. 2 that the PS approach leads to 
the multiplicity function rj(L)dL/L with: 

/ x dim/ [2 po „2 /9 dim/ , . 

^=^dl^ = V^ e ' dinL ^ 
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Fig. 9. The galaxy luminosity function r)(L) dL/L. The graph 
shows the quantity ij(L) as a function of the magnitude in the 
B-band for the non-linear scaling approach (solid line) and the 
PS prescription (dashed line) . The data points are observations 
from Loveday et al.(1992) (circles), Ellis et al.(1996) (filled 
squares) and Zucca et al.(1997) (triangles). 

and v given as a function of M by (2), the mass M being 
in turn related to L by the model of star formation. As we 
could see in Fig. 3, the typical masses of our galactic halos 
are M ~ 2 10 12 M & . As a consequence, since the cutoff of 
the multiplicity function given by the PS approach occurs 
for v ~ 1, which should correspond roughly to the cutoff 
of the observed luminosity function, the number density 
of galaxies is of the order of: 

V (L) ~ g e~ 1 ' 2 ~ 0.03 Mpc" 3 (31) 

which is much larger than the value 77 (L») ~ 3 10 _3 Mpc~ 3 
given by observations. Hence the luminosity function im- 
plied by the PS approach is much larger than the observed 
values at luminosities smaller than its cutoff, as we can see 
in Fig. 9. A galaxy similar to the Milky Way, with a mass 
M ~ 2.6 10 12 M©, a density contrast A ~ 2200 and a B- 
band magnitude Mb — —20.4, corresponds to a linear pa- 
rameter v ~ 1.6. Hence it is already at the cutoff of the PS 
multiplicity function while the cutoff of the observed lumi- 
nosity function corresponds rather to Mb — —21. Thus, 
the luminosity function given by the PS approach falls 
down at luminosities smaller than what is observed, as we 
can see in Fig. 9. Moreover, in the case where P(k) cx k n 
(on galactic scales n ~ —2) we have: 

a(M) cx M- (n+3)/6 (32) 

This allows us to get (App.G) the shape of the luminosity 
function implied by the PS approach in the three regimes 
we have already considered rj(L) cx L~°- b e~ L / Ls , but with 
a rather small value for L s . Thus, as seen in Fig. 9, the 
luminosity function we get in this way has a shape some- 
what similar to the observations (a power-law with an 
exponential cutoff), but, as already discussed by VS, its 
normalization is too high and its cutoff is too strong. 



On the other hand, the non-linear scaling approach 
leads to: 

with x given by (8) and the function h(x) by (29). Now, 
the cutoff at L* must correspond to X X ^ , so. 

»7(L,) ~ g xl h{x.) ~ 10- 3 Mpc" 3 (34) 

with a typical mass M ~ 2 10 12 M©. This is close to the 
observed values, hence the normalization of the luminos- 
ity function implied by the non-linear scaling approach is 
consistent with observations, as we can see in Fig. 9. We 
can also look at the shape of the luminosity function in 
the 3 regimes we have already considered. We get (App.G) 
7](L) cx £-°- 46 at the faint end and r)(L) cx L- 1A e~ L / L - 
for the bright galaxies, with a much larger value of 
than with the PS approximation. Hence we see, Fig. 9, 
that the luminosity function we get in this way is quite 
close to observations: it shows an exponential cutoff for 
bright galaxies Mb < — 21 and a power-law behaviour for 
faint galaxies Mb > —17. The existence of a flat plateau 
for — 21 < Mb < —17 is characteristic of our results (for 
both the PS and non-linear scaling approaches) and holds 
also for other power-spectra (e.g. n — —2). This feature is 
quite remarkable, as it did not appear in previous models 
but is in good agreement with observations, which seem to 
show an upturn at Mb > — 18 after a fiat portion, see for 
instance Driver & Phillipps (1996). This sudden change of 
the slope of the luminosity function around M B — —18 in 
our model corresponds to the transition from the curve Ca 
to C v , which define the global properties of galaxies (mass, 
radius,...). 

As we noticed in Sect. 2 the PS mass function can be 
written in the same form as the non-linear scaling one, 
with a different scaling function h(x) which has a stronger 
and earlier cutoff for large x and a higher normalization at 
x ~ 1 (see also VS). The difference between both luminos- 
ity functions we can see in Fig. 9 is a direct consequence of 
the difference between these scaling functions h(x). This 
suggests in turn that the non-linear scaling function de- 
scribes the actual outcome of gravitational processes more 
accurately than the scaling function derived in the PS ap- 
proach, and hence than the usual PS prescription. Indeed, 
the relation x — L is strongly constrained by the Tully- 
Fisher relation, which thus enables one to derive strong 
constraint on the scaling function h(x). Note however that 
there is an additional factor p /M in the luminosity func- 
tions (33), so that the power-law or the normalization of 
h(x) cannot be constrained without a specific model for 
galaxies, like ours, which gives the value of M attached to 
the parameter x. Nevertheless, the exponential cutoff of 
h(x) is quite strongly constrained by the one of the lumi- 
nosity function. We can note that Kauffmann et al.(1998) 
also found that a Press-Schechter approach predicts too 
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many intermediate galaxies as compared to the results of 
N-body simulations. Although their model is significantly 
different from ours this effect agrees with our analysis of 
the mass functions (see Sect. 2. 3 and VS). 

4.1.8. = 1, n = -2 and n = -1 

Similarly to the case of a CDM power-spectrum we can 
also consider the cases of a power-law P(k) oc k n with 
n = — 2 and n = — 1. The case n = — 2 gives results very 
close to those obtained for a CDM power-spectrum. This is 
natural since the latter has a local slope n ~ — 2 on galac- 
tic scales. On the contrary, the case n = — 1 produces a 
very strong exponential cutoff which leads to a luminosity 
function quite far from the observed one. Indeed, we now 
have a 2 oc M~ 2 / 3 instead of a 2 oc M -1 / 3 (for n = -2) 
which means that bright massive galaxies are much more 
rare relative to small ones as compared to the previous 
case n = — 2, since density fluctuations decrease faster 
with larger mass. Thus a power-spectrum index n = — 1 
at galactic scales seems to be incompatible with the ob- 
served luminosity function, at least within the framework 
of our model. Moreover, as most galaxy characteristics 
are strongly constrained by observations (Tully-Fishcr re- 
lation, gas/star mass ratio, lower limit for galactic masses 
and radii, ...) it is very likely that no reasonable model 
would reconcile the observed galaxy luminosity function 
with n = —1. 

4.2. O = 0.3 , A = 

In the case O = 0.3 and A = 0, we choose Ob = 0.03 and 
£78 = 0.77. The physical properties of the galaxies are close 
to those in the = 1 universe, the analogs of Fig. 3, Fig. 4 
and Fig. 5 show the same behaviour as in the previous 
case. For instance, a galaxy similar to the Milky Way cor- 
responds to R ~ 97 kpc, M ~ 1.6 10 12 M©, A ~ 14000, 
T ~ 3.3 10 6 K and L ~ 3 10 10 L©. Note that we obtain a 
smaller mass and radius, as compared to the case = 1, 
as we explain in App.F. 

We consider the case of a CDM power-spectrum. 
Fig. 10 shows the linear correlation function a(R) 2 and 
the evolved non-linear correlation function at red- 

shift z = 0. 

We see that the evolved non-linear correlation function 
is still reasonably close to a power-law £(R) oc i? -7 , with 
7 = 1.8, in the range 30 kpc to 1 Mpc. The scaling function 
h(x) we need to obtain the galaxy luminosity function in 
the non-linear scaling approach is not available from the 
current numerical simulations but as we argued in VS it is 
expected to be similar to the scaling function obtained in a 
critical universe with the same power-spectrum. Hence we 
adopt the function h{x) we used for a critical universe, see 
(29). The galaxy luminosity function we get in this way is 
shown in Fig. 11 as a function of the B-band magnitude. 




O =O.3, z = 



-2 - 1 O 1 

log(R) [ Mpc ] 

Fig. 10. The evolved correlation function (,(R, z) at z — 
(solid line), and the linear extrapolation a(R) 2 (dashed line), in 
the case Qo = 0.3 , A = and P(k) is a CDM power-spectrum. 
The dot-dashed line is the power-law (R/R )~ J with 7 = 1.8 




Fig. 11. The galaxy luminosity function r)(L)dL/ L. The graph 
shows the quantity rj(L) as a function of the magnitude in the 
B band for the non-linear scaling approach (solid line) and the 
PS prescription (dashed line). The data points are as in Fig. 9. 

As we can see from Tab.l in App.F, when we shift Oo 
from 1 to 0.3 we should decrease Ob to keep the radius 
of the Milky Way larger than 60 kpc. In fact, since for a 
critical universe we had a large radius R ~ 140 kpc, we 
can keep Ob roughly constant. Thus, we choose Ob = 0.03 
which leads to R ~ 97 kpc. From Tab.l it implies that the 
mass of galaxies declines roughly in the same proportion 
as the radius, and indeed we now have M ~ 1.6 10 12 M© 
for the Milky Way. The number density of galaxies given 
by the PS approach is now: 

77(L) ~ e- 1 / 2 ~ 0.01 Mpc" 3 (35) 

which is closer to observations than it was in the case 
O = 1, as we can see in Fig.ll. However, since the rela- 
tion temperature - luminosity did not change (because it 
is constrained by the observed Tully-Fisher relation), the 
cutoff entailed by the PS approach is still too strong. As 
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was the case for a critical universe, the non-linear scaling 
approach can produce a smoother cutoff at higher lumi- 
nosities. The analysis made in Sect. 4. f .7 and App.F for the 
three regimes of galaxies still holds, hence we recover the 
same slopes. However, in such a low density universe the 
curve C v is lower relative to Ca (as compared to a critical 
universe) on the analog of Fig. 3, hence even for relatively 
faint galaxies Mb < —16 we still are in the regime 2), 
thus the slope of the luminosity function is very small, for 
both prescriptions. In fact, as we can see in Fig. 11, the 
luminosity function is flat between Mb = — 16 down to 
Mb = —20, which is quite remarkable. Overall, the agree- 
ment with observations is still very good. 

5. Time-evolution 

The method we presented in the previous sections can 
obviously be used to derive the galaxy properties at any 
time, which allows us to get the evolution with redshift of 
the physical characteristics of galaxies and of their mass 
function or luminosity function. 

5.1. n = 1 

5.1.1. Dark matter properties 

The relation temperature-density which corresponds to 
Fig. 2 is shown in Fig. 12 at the redshifts z — 2, z = 
and z = —0.5. Its behaviour is similar to the one we got 
at z = 0. As we noticed in Sect. 3. 2 the curves C v , hence 
C, depend on z. This changes the shape of C, which will 
modify the shape of the luminosity function rj{L) as we 
shall see. 
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Fig. 12. Temperature-density diagram in the case Q = 1 at 
the redshifts 2 = 2, z = and z = —0.5 from top down to 
bottom. The critical density at z = is noted p c o- As in Fig. 3, 
the short-dashed lines correspond to halos of fixed mass, on 
the left, and halos of fixed radius, on the right. Their position 
does not vary with the redshift. 



We can notice in Fig. 12 that as the redshift increases 
halos get smaller and less massive. However, since Ca does 
not change with z some halos have fixed physical charac- 
teristics below a certain redshift, as long as they stay on 
Ca (in the approximation of stable clustering, which holds 
for £ > 200). We may consider the evolution of halos in 
the extreme regimes 1) and 3): 

1) Faint galaxies located on C v : A = constant. 

At fixed temperature T we get M cx (1 + z)~ 3 / 2 and 
R oc (1 + z)~ 3,/2 . These halos become smaller and less 
massive, and we have: 

^cx (l + z) (1 -™ )/4 and x oc (1 + z) 3 ~ 37/2 (36) 

where we assumed stable clustering. Thus, as long as the 
exponential cutoff plays no role (y < 1 or x < a;*), 
the temperature function given by the PS approach ver- 
ifies n(T) cx (1 + z)( 7 "™)/ 4 , that is r?(T) cx (1 + z) 2 - 25 
if n = —2, while the non- linear scaling approach gives 
r)(T) cx (1 + z) 1 - 6 if 7 = 1.8 and oj = 0.3. Hence the co- 
moving number density of halos at these small tempera- 
tures increases with z in both prescriptions, but somewhat 
more strongly within the framework of the PS approach. 
However, since v and x increase the exponential cutoff will 
eventually lead to a decrease of the comoving number den- 
sity. 

3) Bright galaxies, located on Ca : R= constant 
The physical properties of these halos remain constant 
with the redshift while their density contrast decrease as 
(1 + A) cx (1 + z)~ 3 . This holds as long as A > 177. Their 
parameters v and x (if the clustering is stable) remain 
constant too so the temperature function does not evolve. 
This implies that the changes of the comoving luminos- 
ity function in the range corresponding to these halos will 
only be due to the variation of their luminosities (pure 
luminosity evolution). 

As the redshift increases, the regime of galaxies corre- 
sponding to a fixed temperature changes: 3) — > 2) — > 1). 
Finally, for z > 5 all halos belong to the regime 1). 

5.1.2. Star formation history 

Fig. 13 shows the star formation rate dM s /dt for various 
redshifts as a function of V c . Small galaxies, with a low 
circular velocity, have a small star formation rate because 
their star formation time-scale To is very long (r 3> t) 
since their density and temperature are small. On the 
other hand, very large and luminous galaxies also have 
a small star formation rate in the present universe be- 
cause their star formation time-scale is short as compared 
to their age (to «C t) , so they have already consumed most 
of their gas. In the past they had a higher star formation 
rate since their gas content was larger. 
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log(V„) [ km/s ] 

Fig. 13. Star formation rate in M /year as a func- 
tion of the circular velocity V c in km/s, for the redshifts 
z — 2, z — and z = —0.5 from top down to bottom. 

More precisely, we get for the regimes 1) and 3): 
1) Faint galaxies located on C v : 



We have M a 



Mi 



b oc 



V 3 and t oc V 2 . Since 



dM s /dt ~ Mg/T we obtain dM s /dt oc V" 5 . 

3) Bright galaxies, located on C\ : 

For galaxies which have already consumed most of 
their gas content we have M g oc V 2 cxp(— V/Vo) and 
r oc V- 1 . Hence dM s /dt oc V 3 exp(-V/Vb). 

Fig. 14 shows the time-evolution of the star formation 
rate dM s /dt (upper panel) and of the metallicities Z c , Z h 
and Z s (lower panel) of a galaxy similar to the Milky 
Way. The slight decrease with time of the star formation 
rate after 5 10 9 years is due to the smaller gas content 
of the Galaxy, as its gas is gradually converted into stars. 
This variation by about a factor 2 is consistent with ob- 
servations, which imply in fact that the star formation 
rate of the Milky Way did not vary by much more than 
this amount. Thus, another prescription with a stronger 
evolution would increase the global comoving star forma- 
tion rates or the luminosity functions we shall obtain at 
z <~ 1 — 2, which would improve the agreement with ob- 
servations, but it would not satisfy the mild evolution ob- 
served for the Milky Way. One could build for instance 
such a scenario by using for the star formation time-scale 
r c a different prescription: r c oc l/p g (where p g is the gas 
density in the core) instead of t c oc p~ x l 2 . Such a param- 
eterization, quite plausible, leads naturally to a stronger 
dependence on time, through the density, and to the ef- 
fects we described above, as we checked numerically. To 
keep things simple we shall keep our original prescription, 
which satisfies the constraints given by the Milky Way. 
The slight increase of the star formation rate shortly after 
the time p t m is due to the fact that some gas is still falling 
onto the galaxy inner parts from the halo. 
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Fig. 14. Upper figure: evolution with time of the star forma- 
tion rate dM s /dt of the Milky Way. Lower figure: evolution of 
the metallicities Z c , Zh and Z 3 of the Milky Way. The vertical 
dashed- line represents the time p t m . The marks on the upper 
left side of both figures show the redshift. 

The temporary increase with redshift of the star for- 
mation rate at small times t < p t m (on the left of the 
vertical dashed-line) is due to the decrease of the relevant 
time-scales r c and r<2 at high redshift: r c <~ oc p^ 1 / 2 oc 
(1 + z)~ 3 / 2 . Note however that during this epoch the mass 
which will later form the Milky Way is distributed over 
several smaller halos. The metallicities (meant as [O/H] 
rather than [Fe/H] since we do not include SNIs) increase 
steadily with time, in all three components (M gc , M g h and 
M s ). At the time p t m , the metallicities Z c and Z s were 
smaller than their present values by a factor A[Z] = —0.5, 
which agrees with the difference in metallicity between the 
oldest and youngest stars in the disk. Stars created before 
this date with a high metallicity —1 < [Z] < —0.5 are lo- 
cated within the bulge, while even older and less metallic 
stars would be located in the halo. The very steep increase 
in metallicity at the begining implies that there are very 
few stars of very low metallicity [Z]< —2. 

Fig. 15 shows the cumulative stellar metallicity distri- 
bution P(< Z c ) for the Milky Way: 



P(< Z) 



M lo (t z ) 



with Z c (tz) = Z 



(37) 



Mi o (t ) 

Thus, at time p t m , which corresponds to the formation 
of the disk, the metallicity of the gas, and of the stars 
being formed, was a factor 0.3 smaller than the metallic- 
ity of stars created today. Moreover, 25% of present long- 
lived stars were formed by this date, which is roughly the 
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Fig. 15. The cumulative stellar metallicity distribution for the 
Milky Way. The vertical dashed-line represents the metallicity 
of the newly formed stars at the time p t m , when the object 
switches from the merging regime to the isolated galaxy regime. 

stellar mass of the bulge. Thus, in our model bulges of 
disk galaxies contain an important proportion of old stars 
(with Z ~ 0.3Z Q for the Milky Way) which formed be- 
fore the disk along processes similar to those of elliptical 
galaxies. This agrees with the interpretations of Ortolani 
et al.(1995) (who found bulge globular clusters as old as 
halo globular clusters) and of Jablonka et al.(1996). Note 
however that some amount of star formation may have 
kept going on until today. The slow decrease of the metal- 
licity in Fig. 15 with the fraction of the mass of stars also 
shows that we have no G-dwarf problem. 



— —5 





1 


iii 


1 1 


























5 \ 


V 










\ 










\ 








1 ij I L_ 


\ a 
\ 

\ 

i \, 




1 







2 






log(dM,/dt) [ 


M / year ] 







Fig. 16. The distribution function of the star formation rates 
at different redshifts: z = (solid line), z = 1 (dashed line) , 
z = 2 (dot-dashed line) and z — 5 (dotted line) . 

Fig. 16 shows the distribution function of the star for- 
mation rates at different redshifts (i.e. the number of 
galaxies with a given star formation rate per comoving 
Mpc 3 ). It extended to much higher star formation rates 
in the past, at z ~ 1 — 2, than it does in the present uni- 
verse, because massive galaxies formed most of their stars 



at these early epochs (since they have a small star for- 
mation time-scale), when they experienced a very active 
phase, and their star formation rate has steadily declined 
ever since as their gas content became smaller. This cor- 
responds to the history we developed previously in detail 
for a galaxy similar to the Milky Way. Moreover, the lumi- 
nosity (or the galactic mass) of the galaxies characterized 
by the highest star formation rate was larger in the past 
(at z ~ 1) than it is now. These results agree with the 
redshift evolution observed by Cowie et al.(1996, 1997). 

5.1.3. Luminosity evolution 

We can see in Fig. 17 that when the redshift increases 
short-lived stars become more important because these 
halos arc younger, have more gas and a higher star for- 
mation rate. Hence the global mass-luminosity ratio gets 
smaller and, at fixed V c , massive halos which are on the 
cooling curve C\ have a larger luminosity since their mass 
remains constant (as long as they remain on C\). Thus, 
the slope of the temperature-luminosity relation at high V c 
gets stronger. As a consequence, the knee of the luminos- 
ity function should move toward larger luminosities in the 
past, and fainter luminosities in the future, since we no- 
ticed above that the temperature function in this region 
docs not evolve with z. On the other hand, at small V c 
the temperature-luminosity relation keeps the same slope 
since the analysis developed for z = is still valid. 

However, the normalization decreases at higher red- 
shifts, although the mass/luminosity ratio is smaller, be- 
cause the mass of these halos decreases. Indeed, at higher 
redshifts the time r*, which is the life-time of stars at 
the boundary between our two classes of stars (short-lived 
which are recycled and long-lived ones), scales as t* oc 
iff oc (1 + z) -3 / 2 , see App.D. The mass/luminosity ratio 
of the global stellar population varies more slowly than 
r» because the IMF contains fewer massive stars: M s /L cx 

T l-(x-l)/(,-l) K (1 + z) -3/2+3/2 where x and „ 

are stellar parameters (IMF, mass- luminosity relation, see 
App.D). Since for faint galaxies, located on C v , we have for 
a fixed temperature M s oc M oc (1 + z)~ 3 / 2 , because the 
global star formation time-scale r follows the decrease of 
the galactic age, the luminosity decreases slowly at high 
redshifts, as we can check in Fig. 17. As galaxies leave the 
regime 3) to enter the regimes 2) (even when they remain 
on Ca, the change is that now their stellar content is small: 
M s < M g and M g ~ M b ) and finally 1), they satisfy the 
scaling V c — L described in Sect. 4. 1.3 and App.G for the 
regime 1) at z = 0. This explains why the slope of the 
Tully-Fisher relation we get remains constant with z, and 
extends up to the bright galaxies at z > 2 (the high lumi- 
nosity bend we have for z < disappears). The regime 1), 
which only corresponded to the smallest galaxies at z = 
is now valid for nearly all galaxies at z > 2. 
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Fig. 17. Upper figure: ratio of the total mass of stars M s over 
the galaxy luminosity L, as a function of the circular velocity 
V c , for the redshifts z = 2, 2 = and z = —0.5. Lower figure: 
luminosity L of the galaxy as a function of V c . The dot-dashed 
curves show the luminosity due to massive stars (L s h M a h). At 
large V c a higher z corresponds to larger L and larger L s h M s h ■ 

5.1.4. Luminosity function 

Fig. 18 shows the evolution with redshift of the comoving 
galaxy luminosity function in the case of a CDM power- 
spectrum. We can see that the knee of the luminosity 
function moves toward larger luminosities in the past, un- 
til z = 1, and fainter luminosities in the future, while 
the comoving number density for faint galaxies increases 
until z < 2. This is consistent with what we expected 
from the above analysis. Indeed, for faint galaxies we saw 
previously that r/(T) oc (1 + z) 2,25 for the PS approach 
and rj(T) oc (1 + z) 1,5 for the non-linear scaling prescrip- 
tion. Since the Tully-Fisher relation does not change very 
much in this range, the luminosity function r/(L) follows 
the same behaviour. 

On the contrary, for bright galaxies we explained above 
that while the temperature function rj(T) does not evolve, 
the slope of the Tully-Fisher relation gets larger, so that 
a given temperature corresponds to a larger luminosity in 
the past, hence the luminosity function cutoff moves to- 
ward larger luminosities. We can also notice that the faint- 
end slope of the luminosity function gets higher in the 




Fig. 18. The comoving galaxy luminosity function r)(L) dL/L 
at the redshifts 5, 2, 1, and —0.5, for the non-linear scaling 
approach (solid line) and the PS prescription (dashed line). 
The data points are observational results at z = as in Fig. 9. 

past. As we said previously, this is related to the change 
of the curve C with z. Indeed, these faint galaxies leave 
the regime 2) to enter 1), hence the slope of the luminos- 
ity function increases from rj(L) oc L~ 2 to r](L) oc L -0 ' 5 
for the PS approach, and from r](L) cx _L~ 24 to rj(L) oc 
£-0.46 £ or £ ne non-linear scaling prescription (note that 
although this increase is qualitatively correct, as we can 
see in Fig. 18, the slopes of the luminosity function we ob- 
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tained in this way are not very accurate because we did 
not consider the variation of the ratio M s /L). This ef- 
fect could explain the steepening of the faint-end slope of 
the luminosity function which is observed in the past at 
z = 0.5 or z = 1. Moreover, the slope of the luminosity 
function after the cutoff gets smaller, especially for z > 2. 
Indeed, we still have (with n ~ —2): 

V (L) OC M^' 6 e -(M/M,)<"+s>/3 or ^ K MWs e -M/M, 

but as galaxies leave the regime 3), where M g <C M s ~ 
Mf,, to enter the regime 2), where M s -C M 9 ~ M(,, as 
their gas content increases, the mass-luminosity relation 
becomes L oc M 5 / 2 instead of L oc M. Hence we get: 

rj(L) oc Z/(™ -1 )/ 15 e -(V^) 2< " +5)/15 p S appr o a ch 

?y(L) cx £ 2w =/ 5 e ~( L / L *) 1 non-linear scaling approach 

With n = — 2 and w s = —3/2 we obtain: 

V (L) oc i-°' 2 e -W L ')" 4 or r?(L) oc iT - 6 e -(W" 

Thus, in both cases the exponential cutoff becomes less 
sharp in the past (note that the cutoff characteristics 
or L, are not the same for the PS and non-linear scal- 
ing approaches, as we noticed earlier, in Sect. 4.1.1. for 
instance) . 

We can see that the galaxy comoving number density 
evolves much faster for the non-linear scaling prescription 
than for the PS approach. Indeed, it leads to a luminos- 
ity function which is much higher than the PS one at the 
bright end for z < 2 and it suddenly decreases at z > 2 to 
superpose onto the PS prediction at z ~ 4 and then gives 
even fewer galaxies than this latter prescription. The slow 
evolution for the PS case is due to the fact that v, hence 
r)(T), is constant in the regime 3), and only increases as 
(1 + z)^ 1- ™)/ 4 once the halo enters the domain 1). More- 
over, when v is not too large, in this regime 1) the pref- 
actor grows as (1 + z)^ - ™-*/ 4 which balances the decrease 
of the exponential term. On the contrary, in the case of 
the non-linear scaling prescription x increases slowly until 
£ ~ 177, since the clustering is not exactly stable. Then, 
when £ < 177 the correlation function £ decreases sud- 
denly very strongly, and x rises sharply. This produces a 
sharp decrease of the comoving number density, as the ex- 
ponential cutoff becomes very important. We must note 
that for £ <C 177 the function h(x) should change so that 
for £ <C 1 it becomes equal to the result obtained in the 
quasi-gaussian regime (see Colombi et al.1997). However, 
the exponential cutoff of the the quasi-gaussian function is 
stronger than for the non-linear case (Colombi et al.1997), 
so the predicted decrease of the luminosity function would 
be somewhat stronger. Thus, it is important to note that 
up to z ~ 3, which is already a rather large redshift, the 
non-linear scaling approach leads to a much higher lumi- 
nosity function than the PS prescription at the bright end. 



In this regime, £ is still of the order of ~ 10 at the onset 
of the exponential fall-of. At such values of £ it has been 
seen (Bouchet et al.1991) that h(x) is still given to a good 
approximation (10%) by its non-linear form. 

5.1.5. Average comoving stellar properties 

The evolution with redshift of the stellar density param- 
eter ft s (z)/£lb, see (Cll), is displayed in Fig. 19 (upper 
panel). Its decrease in the past is due to two effects. First, 
we only consider halos with a temperature T > 10 4 K 
(cutoff due to inefficient cooling), which leads to higher 
parameters v and x at high redshifts, see (36), hence to 
a smaller mass fraction - there is less mass contained in 
deep potential wells in the past. Second, at larger redshifts 
galaxies have a higher gas/star mass ratio (smaller t/ro), 
hence most of the mass is in the form of gas. Note that 
in the present universe, since only galaxies more luminous 
than the Milky Way have an appreciable stellar content 
(our Galaxy is just between both regimes as M g /M s ~ 1) 
and they form a small part of the total mass (they already 
are in the exponential cutoff of the luminosity function) 
the ratio Cl s (z — 0)/fi{, is still small. 




z 



Fig. 19. Upper figure: evolution with redshift of the stellar 
density parameter Q s (z)/flb- Lower figure: comoving star for- 
mation rate dp a /dt for < z < 5. The solid lines correspond 
to the case fi = 1 and the dashed lines to Qo — 0.3. In both 
cases we only display the non-linear scaling prescription. The 
data points are taken from Madau et al.(1996) (squares), Lilly 
et al.(1996) (disks), Gallego et al.(1995) (triangle) and Cowie 
et al.(1995) (crosses). Note that the points at high redshift 
(z > 2.5) are only lower limits. 
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We can see in the figure that two thirds of the present 
mass in stars formed recently at z < 1. It is also clear 
that we do not encounter the usual overcooling problem - 
all the gas cools and is converted into stars within small 
objects at high redshifts because cooling is very efficient 
(high densities) - within this framework. As we explain 
in App.C.2, this is due to the redshift dependence of our 
star formation time-scale To, which ensures that although 
the gas may cool at high redshift it cannot be immedi- 
ately converted into stars (to does not decrease faster than 
the age of the universe with redshift). Hence there is still 
plenty of gas available in the present universe (£l s <C f2&) 
which is not necessary cold as it has been reheated by su- 
pernovae, stellar winds and by the energy released by halo 
mcrgings and collapse. 

The comoving star formation rate, see (C13), is shown 
in Fig. 19 (lower figure). It first increases with redshift un- 
til z <~ 1 because the star formation time-scale To de- 
creases as To oc (1 + z)~ 3//2 , as long as T > Tq. However, 
at high redshifts z > 1 the comoving star formation rate 
gets smaller because the mass contained in deep poten- 
tial wells starts to decrease and To gets larger because of 
the factor (1 + Tq/T). Nevertheless, its variation over the 
whole range < z < 5 is rather small. 



5.1.6. Galaxy counts 

Finally, the evolution of the luminosity function allows us 
to calculate the galaxy number counts as a function of the 
apparent magnitude tub and the redshift z. The absolute 
B band magnitude Mb is related to ms by: 



m B = M B + 5 log 



ri(z) 



lOpc 



K(z) + E(z) 



(38) 



where rt(z) is the luminosity distance to redshift z, K{z) 
is the usual K-correction and E(z) is the evolution correc- 
tion. These last two terms vary with the stellar and mor- 
phological properties of galaxies, and should be evaluated 
from our galaxy evolution model, to get a self-consistent 
result. However, since this would require a detailed de- 
scription of the stellar properties of galaxies, including 
their spectra and colors, which we plan to tackle in a future 
paper, we shall simply take in this article K(z)+E(z) = 0. 
Indeed, as can be seen in King & Ellis (1985) for instance, 
while K(z) is positive because of the shape of the spec- 
trum, E{z) is negative because galaxies where bluer in 
the past, and both terms cancel roughly. Of course, the 
net result varies with the galaxy type, and the evolution 
model used to get E(z), and for faint magnitudes where 
the contribution to galaxy counts extends to high redshifts 
[z <~ 1) we may have an error of one magnitude. 

However, this approximation should give a fair idea 
of the implications of our model on the galaxy num- 
ber counts, and it does not influence the comparison be- 
tween the PS and non-linear scaling approaches. We note 
J\f(rriB, z)drriBdz the number of galaxies per square degree 
with apparent magnitude mj to wb I dms and redshift 
z to z + dz: 



Af(rriB,z) dms dz 



1 180/ 



${M B ,z) 



dV 

duidz 



dm B dz(39) 



where f (Mg, z)dMs = i](L)dL/L is the comoving lumi- 
nosity function at the redshift z. Finally, the number of 
galaxies N(mB)drriB per square degree is the integral over 
z of A/"(ms, z)dvriBdz: 



Fig. 20. The B-band comoving luminosity density. The solid 
line corresponds to the case f2 = 1 and the dashed line to 
fio = 0.3. In both cases we only display the non-linear scaling 
prescription. 



N(m B ) 



= (—)' 
V180/ 



*{M B ,z) 



Thus, the PS approach leads to: 



dV 
dtodz 



dz 



(40) 



The evolution with redshift of the B-band comoving 
luminosity density is shown in Fig. 20. The luminosity den- 
sity decreases in the past because the mass fraction within 
galaxies diminishes and their stellar content is lower. How- 
ever, the luminosity density remains high until z ~ 3, 
when it is still equal to one half of its present value. 

On all these figures, we can see that the evolution is 
slower for a low-density universe, as is well-known. 



N{m B ,z) 



fJLY 



dV 



2 po 

7T M 



z P0 -v 2 /2 

e ' 



8v 



V180/ dcodz V 7T M 8Mb 
while the non-linear scaling prescription gives: 



Af(m B ,z) 



/_n_\ 2 dV po_ 8x 
1 180/ cLudz M X [X> 8Mb 



(41) 



(42) 



Fig. 21 shows the B counts for the PS and non-linear 
scaling prescriptions, as well as for a no-evolution model 
such that the comoving luminosity function does not 
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Fig. 21. The B band differential number counts. The graph 
shows N(m,B) (number of galaxies per square degree per appar- 
ent magnitude) for the non-linear scaling approach (solid line), 
the PS prescription (small dashed line), and a no-evolution 
model (dot-dashed line) such that the comoving luminosity 
function does not change with z. The data points are taken 
from Lilly et al.(1991) (squares) and Metcalfe et al.(1991) (tri- 
angles) . 

change with time and is equal to the one given by the non- 
linear scaling prescription at z = 0, which allows to distin- 
guish the effects of evolution from those due to geometry 
(comoving volume clement, luminosity distance). Our re- 
sults agree reasonably well with observations although the 
number counts given by the non-linear scaling model are 
somewhat too small at very faint magnitudes ms > 26. 
The counts given by the non-linear scaling approach are 
lower than those obtained by the non-evolving model be- 
cause the number of L» galaxies we get with our model 
decreases when we look in the past (see Fig. 18) and for a 
fixed apparent magnitude the counts are mainly sensitive 
to the evolution of the number density of bright galaxies, 
since we do not see any longer the faintest galaxies. This 
discrepancy with the observational data for tub > 26 is 
a well-known problem for non-evolving models. In fact, 
Cole et al.(1992) showed that natural models where the 
comoving luminosity function evolves in a homogeneous 
way (evolution of the normalization or of the cutoff) so 
as to match the B counts will contradict the data in the 
K band and the observed redshift distribution of galaxies. 
Hence they concluded that a new population of rapidly 
evolving blue galaxies fainter than L s is necessary to fit 
all the data. Such an effect, even if real, cannot be given by 
our simple parameterization of star formation. This point 
deserves a detailed study that will be done elsewhere. 

Fig. 22 shows the redshift distribution J\f(rriB,z) of 
galaxies selected at an apparent magnitude mj = 24.15, 
for the PS and non-linear scaling prescriptions, and for 
the no-evolution model. Of course we recover the same 
features as for the integrated number counts: the PS pre- 
scription (because of its high normalisation) and the non- 
evolving model give more galaxies than the non-linear 



Fig. 22. Redshift distribution of galaxies with apparent B 
band magnitude m_g = 24.15. The graph shows 7V(24.15, z) for 
the non-linear scaling approach (solid line), the PS prescription 
(small dashed line), and the no-evolution model (dot-dashed 
line). The data points (histogram) are from Cowie et al.(1996) 
for galaxies in the range 22.9 < ms < 24.4. 

scaling approach. For these last two cases we find a peak 
at z ~ 0.3 which is consistent with observations (Cowie et 
al.1996; Colless et al.1993; Colless et al.1990; Broadhurst 
et al.1988) for this apparent magnitude. Note also the fast 
evolution of the PS result. 

5.2. n = 0.3 , A = 

In the case flo = 0.3 and A = the analysis developed 
for a critical universe still holds. Fig. 23 shows the evolu- 
tion of the comoving galaxy luminosity function with the 
redshift in this case for a CDM power-spectrum. As we 
can see, the evolution is qualitatively similar to what we 
could see in Fig. 18 but much slower. This is a well-known 
property of low-density universes: the evolution of gravi- 
tational clustering is slower than for a critical universe. 

However, the quantitative difference with the case = 
1 for the bright end of the luminosity function is quite 
dramatic, and thus appears to be extremely sensitive on 
the cosmological parameter ri . 

Fig. 24 shows the B counts for the PS and non-linear 
scaling prescriptions, as well as for the no-evolution model. 
As was the case for O = 1, the slope of the counts gets 
smaller at faint magnitudes raj > 25. 

Fig. 25 shows the redshift distribution of galaxies se- 
lected at the apparent magnitude raj = 24.15. Once again 
we can see that the evolution is slower than for £1 = 1 : the 
peak moved to higher redshifts z ~ 0.4. 

6. Quasar number density 

As quasars are among the first objects to form they pro- 
vide an interesting probe of the Universe at high redshifts 
z ~ 2 — 5. Moreover, they constrain the models of grav- 
itational clustering to form objects of large mass at high 
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Fig. 23. The comoving galaxy luminosity function rj(L) dL/L 
at the redshifts 5, 2, 1, and —0.5, for the non-linear scaling 
approach (solid line) and the PS prescription (dashed line). 
The data points are observational results at z — as in Fig.ll. 

redshifts. We shall now see whether our formalism can 
satisfy this requirement. We first assume that the quasar 
mass Mq is proportional to the mass of gas M gc available 
in the inner parts of galaxies: 

M Q = F M gc (43) 

For galaxies which have not already transformed most of 
their gas content into stars we can use (D5). With (B5) 
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Fig. 24. The B band differential number counts. The graph 
shows AT(ms) (number of galaxies per square degree per appar- 
ent magnitude) for the non-linear scaling approach (solid line), 
the PS prescription (small dashed line), and a no-evolution 
model (dot-dashed line) such that the comoving luminosity 
function does not change with z. The data points are as in 
Fig.21. 




2 



z 



Fig. 25. Redshift distribution of galaxies with apparent B 
band magnitude m B = 24.15. The graph shows 7V(24.15, z) for 
the non-linear scaling approach (solid line), the PS prescription 
(small dashed line), and the no-evolution model (dot-dashed 
line). The data points are from Cowie et al.(1996) for galaxies 
in the range 22.9 < m B < 24.4. 

this means that the quasar mass is proportional to the 
stellar mass Mq ~ F M s . We write the bolometric lumi- 
nosity Lq of the quasar as: 

L V = (44) 

where e < 0.1 is the quasar radiative efficiency (fraction of 
central rest mass energy converted into radiation) and £q 
is the quasar lifetime. If quasars radiate at most at the Ed- 
dington limit at which radiation pressure on free electrons 
balances gravity (Efstathiou & Rees 1988; Nusser & Silk 
1993) one has t Q8 > 4.4 e where t QS = t Q /10 s yr. Finally 
we use a bolometric correction factor Lb/L^i = 0.16. 
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Since we only consider here very massive and rare quasars 
we write the quasar comoving number density as: 



Vq(Mq) 



dM Q 



M, 



Q 



I [v(M,z) 
d 

ft Q —r](M,z) , 
y dt x M 



v(M,z[t(z) 
dM 



t Q ])] 



dM 
~M 

(45) 



in a fashion similar to Efstathiou & Rees (1988) and 
Nusser & Silk (1993). Here r\{M,z) is the galaxy co- 
moving multiplicity function and we assumed a fraction 
/ < 1 of galaxies actually contains a quasar. Thus we 
only have two parameters: (e F/tQs) (Qb/Qo) which sets 
the quasar mass, from (43) and (44), and (/ igs) which 
enters the multiplicity function in (45). Thus a first esti- 
mate of the comoving number density of quasars brighter 
than M B = -26.7 (L B > 8 10 12 L ) is: 



dz 

N Q1 (> L B ,Z) =ft Q -^ fo N (> M rmn,z) 



(46) 



where N(> M m i n , z) is the number of galactic halos more 
massive than M m i n at the redshift z. The minimum galaxy 
mass M min is given by (44). However, the mass of gas 
available in very massive galaxies does not increase lin- 
early with the dark matter mass of the parent halo be- 
cause in these very dense galaxies star formation was very 
efficient so that they have already consumed most of their 
initial gas content, as can be seen from (B6) and (B7). 
As a consequence, the mass of gas which may power the 
central black hole reaches a finite maximum for galaxies 
somewhat more massive than the transition between C v 
and C\. We can note that the maximum quasar mass (or 
luminosity) obtained in this way decreases as time goes on, 
along with the decline in the mass of gas which is progres- 
sively turned into stars. This very simple effect will obvi- 
ously have important consequences, at variance with the 
predictions of models like (46) where the typical quasar 
mass keeps increasing with time as larger scales become 
non-linear. Thus we write a second estimate of the quasar 
comoving number density as: 

N Q2 (> L B ,z) = N Q1 (M mm ,z)~N Q1 (M max (z),z) (47) 

where M max [z) is three times larger than the mass which 
corresponds to the intersection of the curves Ca and C v 
at the redshift z. This takes care of the upper bound 
on the possible quasar mass, as described above. When 
M max (z) ^> M min (which is always the case at high red- 
shifts) the cutoff at M max (z) has no influence and we re- 
cover the number density of quasars given by the previous 
calculation (46). 

6.i. n = i 

In the case = 1 with a CDM power-spectrum, Fig. 26 
shows the comoving number density of bright quasars Nq 
as a function of redshift for both prescriptions: Nqi(z) and 



Nq 2 (z). We use (e F/t QS ) = 2.3 1(T 3 and / t Q8 = 0.018. 
This corresponds for instance to F = 0.01, e = 0.1, tQs = 
4.4 e (Eddington luminosity) and / = 0.04. Then, the min- 
imum dark matter mass for M B < —26.7 is 3.7 10 12 M Q . 
Note that for nearby galaxies Mq <~ 0.006 M s t (Magor- 
rian et al.1998) where M s b is the stellar mass of the bulge 
(assumed to have formed at the same time as the central 
black hole) and that M gc may be smaller than M s by a 
factor 2 for the bright galaxies we consider here, see (D4). 




Fig. 26. The comoving number density of bright quasars 
Nq\(z) and Nq 2 (z). These two estimates superpose at high 
redshift and Nq2(z) shows a sharper cutoff at low redshift. 
The observational points are from Pei (1995). 

As exlained above, both prescriptions Nqi(z) and 
Nq2{z) superpose at high redshift and match the data 
for z > 2. We recover the broad maximum between 
2 < z < 3.5 obtained by observations (see also Hartwick & 
Schade 1990 and Warren et al.1994). At low redshifts, the 
number density of bright quasars given by Nqi (z) declines 
too slowly as compared with observations. On the other 
hand, Nq 2 {z) predicts a very fast decrease so that there 
are no more bright quasars at z < 1. Of course, observa- 
tions show some quasars at z < 1 but this smoother cutoff 
could certainly be obtained with a more detailed model 
which would include some scatter in the mass-luminosity 
relation (and in the properties of galaxies). Moreover, the 
physics of quasars is certainly much more complex than 
the modelisation we used here. Hence we think that these 
results show that our description is consistent with obser- 
vations (a better agreement with the data would require a 
more refined model of quasars themselves to be meaning- 
full) and provides a very natural simple model for quasars 
and galaxies. Note also that our prescription for quasars 
is a straightforward by-product of our model for galax- 
ies: it is its simplest possible extension and we did not 
have to introduce an ad-hoc redshift dependance in the 
mass-luminosity relation to obtain a low redshift decline 
contrary to Haehnelt & Rees (1993). Note that the counts 
at high z are very sensitive to the normalization of the 
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power-spectrum and to the mass of the dark matter halos 
associated to quasars. 

6.2. n = 0.3 , A = 

Fig. 27 shows the comoving quasar number density of 
bright quasars Nq as a function of rcdshift in the case 
O = 0.3, A = for a CDM power-spectrum. We now 
use F = 0.007, e = 0.1, t QS = 4.4 e and / = 0.07. This 
leads to M min = 2.1 10 12 M Q . Thus the mass of the dark 
matter halo which corresponds to a given quasar lumi- 
nosity is lower because the fraction of baryonic matter in 
this universe is larger. As we can see in the figure we ob- 
tain a behaviour similar to the case = 1. However, the 
evolution with redshift of the number density of virializcd 
halos is much slower which implies that the decline of the 
quasar comoving number density at high z is slower than 
for a critical universe (the same effect also appeared for 
the galaxy luminosity function). 
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Fig. 27. The comoving number density of bright quasars 
Nqi(z) and Nq2{z). These two estimates superpose at high 
redshift and Nq2(z) shows a sharper cutoff at low redshift. 
The observational points are from Pei (1995). 



7. Conclusion 

In this article we build a model of galaxy formation and 
evolution based on the hierarchical clustering picture. 

The counts of non-linear objects are done directly 
at the epoch under consideration using the actual non- 
linear density field (instead of its linear extrapolation) by 
methods (Schaeffer 1984; Balian & Schaeffer 1989) that 
have been shown to represent accurately the density field 
obtained in numerical simulations (Bouchct et al.1991; 
Colombi et al.1997). The result, discussed in detail in a 
previous paper (VS), bears some qualitative resemblance 
with the Press-Schechter prescription (which explains why 
the latter appears to be accurate in some cases), but dif- 
fers quantitatively. Indeed, the Press-Schechter approxi- 



mation predicts fewer extreme objects (i.e. with a den- 
sity contrast much larger or smaller than the average fluc- 
tuation at fixed mass, or with a mass much smaller or 
larger than M* at fixed density threshold A ~ 200 where 
cr(M*) = 1) and too many intermediate halos (VS). More- 
over, our approach allows us to go beyond the traditional 
Press-Schechter prescription. Thus we consider mass func- 
tions of objects which are defined by a density contrast 
A(M) which depends on their mass M (see also Valageas 
et al. 1998a). 

Indeed we do not assume that galaxies are defined by 
a constant density threshold resulting from the virializa- 
tion condition, since with such an assumption one would 
count induely virialized clusters as galaxies (Schaeffer & 
Silk 1985). We add to the virialization condition a cooling 
constraint (Silk 1977; Rees & Ostriker 1977) which sepa- 
rates galaxies and clusters. It states that in order to form 
a galaxy a non-linear object must see its baryonic content 
cool and settle in its central part before it gets embedded 
into a larger non-linear halo. More precisely, galaxies are 
the largest patch which satisfies these criteria. Virialized 
objects which do not fulfill this constraint become groups 
or clusters of galaxies. Note moreover that these groups 
are divided into several distinct galaxies which we take 
into account as individual galactic halos (which satisfy the 
cooling condition). Then, we can see that in the present 
universe small and faint galaxies are still defined by the 
virialization condition, which is more restrictive than the 
cooling constraint for these small virial temperatures, and 
are thus the result of many mergers which have kept going 
on until the present epoch. On the contrary, massive L* 
galaxies are given by the cooling constraint, which shows 
that the merging process is no longer effective and their 
evolution is determined by internal physical processes. At 
earlier epochs (higher redshifts), the proportion of objects 
which follow the virialization condition, and are thus the 
result of continuous merging processes, gets larger. The 
£* galaxies reach this regime at z ~ 2, a value which is 
determined solely by the density threshold corresponding 
to virialization (A c ~ 200) and the cooling processes. This 
epoch represents the transition from a regime of continu- 
ous mergers to a regime where galaxies evolve quietly and 
can form a disk. 

To obtain the luminous properties of galaxies we add 
to this evolution of the number of galactic halos a model 
of star formation. We use a star formation time-scale pro- 
portional to the dynamical time (which is also the period 
of revolution of a proto-stellar cloud around the galaxy) . 
In addition supernovae winds are assumed to slow down 
the star formation rate as they can eject or heat some gas 
(White & Rees 1978; Dekel & Silk 1986; White & Frenk 
1991). To this purpose we construct a two-components 
model where the gas is expelled to the outskirts of the 
dark halo with an efficiency that is large when the poten- 
tial well is shallow (which corresponds to faint galaxies) 
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and is not available for star formation until it eventually 
falls back into the baryonic core. This ejected gas may be 
the source of the metal-rich intergalactic gas in clusters, 
which we shall study in a forthcoming article (Valageas 
& Schaeffer 1998). All unavoidable parameters associated 
with star formation processes are adjusted to the well- 
known properties of the Milky Way (mainly its luminosity 
and the gas fraction in the solar neighbourhood), and not 
to the luminosity function to be reproduced. 

Then, our model reproduces the observed Tully-Fisher 
relation, the characteristic luminosities, masses and radii 
of galaxies. Linked with our prescription for the galaxy 
mass function, it provides a prediction for the galaxy lu- 
minosity function which is consistent with observations 
for a critical universe Q = 1 and a low density universe 
£!o = 0.3 , A = with a CDM (or n = —2) initial power- 
spectrum. However, in this latter case the lack of computa- 
tional results to which our non-linear scaling model could 
be compared makes our prediction less certain. On the 
other hand, this very model with the same initial power- 
spectrum and star formation processes shows important 
discrepancies with the observed luminosity function when 
we use a Press-Schechter-like mass function. This discrep- 
ancy cannot be healed without giving away the agreement 
with the above observations, and was well-known in this 
type of approach (inconsistency with the Tully-Fisher re- 
lation, too many small galaxies, incorrect color evolution 
with mass). This, obviously, is no longer a problem in our 
approach, the improvement arising from our better un- 
derstanding of the mass function and our more logical as- 
sumption that baryon concentrations separate when they 
are too large for cooling to occur within a Hubble time at 
formation. 

One interesting new feature which arises from this 
quite classical model of star formation, but less classical 
with respect to the mass function and the definition of 
galaxies, is that massive galaxies are more evolved and 
redder than small ones, as is observed, although our clus- 
tering picture is entirely hierarchical. This was not the 
case for earlier hierarchical models which often used a con- 
stant density threshold to define objects. For this reason, 
it is often stated that hierarchical models do not lead to 
the correct colors of galaxies and have difficulties to pro- 
duce very bright and red galaxies. This statement is shown 
here not to be related to the hierarchical clustering pic- 
ture, but to the modelling of galaxies themselves. Note 
that dust extinction and metallicity effects (not included 
in this article) would also make bright galaxies redder than 
would be obtained otherwise (Kauffmann & Chariot 1998; 
Somerville & Primack 1998). The metallicity produced by 
our star formation process is consistent with the observed 
correlation with luminosity. It is interesting to note that 
for the Milky Way the metallicity of stars formed during 
the early merging regime at z > 2 is \Z\ < — 0.5 which is 



indeed the metallicity cut between disk and halo stars. 

Finally, our analytic model also enables us to derive ex- 
plicit approximate scaling relations between various prop- 
erties of galaxies. This clearly shows the consequences of 
the phenomenological prescriptions adopted for star for- 
mation or feed-back processes for instance, as well as the 
variation of galactic characteristics with the cosmological 
parameters. In addition, we recall here that our model is 
part of a global description of structure formation in the 
universe, which deals with clusters, Lyman-a clouds and 
rcionization processes (detailed in other articles) as well 
as with the galaxies we studied in this article. Thus, our 
goal is to provide a simple realistic analytic model which 
can describe in a unified consistent fashion these various 
phenomena. Despite the simplifications envolved in such 
an approach, this allows one to obtain reasonable quan- 
titative estimates of these physical processes, to draw an 
explicit link between these different aspects of the same 
underlying density field and to get a complementary tool 
to numerical calculations. 

Our approach differs from earlier works by many as- 
pects. White & Frenk (1991) use the Press-Schcchter ap- 
proximation to count overdensities with a fixed density 
threshold, and within each of these halos consider the cool- 
ing condition to determine the fraction of gas that is able 
to cool, arguing the latter will eventually turn into stars. 
Cooling, in this approach, does not modify the number 
of objects, but simply the properties of the latter. Thus 
they do not solve the old problem (Schaeffer & Silk 1985; 
White et al.1987) of predicting too many low-mass galax- 
ies and extremely massive large galaxies (more precisely, 
they only obtain a "halo luminosity function" where the 
most massive objects are not galaxies but groups or clus- 
ters). For similar reasons, in order to avoid the forma- 
tion of too many huge galaxies Kauffmann et al.(1998) 
and Somerville & Primack (1998) impose an ad-hoc cutoff 
(V c < 500 km/s or V c < 400 km/s) to the halos which are 
allowed to cool. This problem is cured in our model by the 
separation of galaxies and galaxy clusters by means of a 
mass-dependent density threshold as implied by our new 
implementation of the cooling constraint. The higher den- 
sity threshold insures the proper falloff at a few 10 12 M Q 
and the slope around 10 6 K of the cooling curve (see Fig. 2) 
provides for the flat slope of the luminosity function. We 
still get a rather steep slope for very small luminosities, 
similar to the one observed (see comments by Driver & 
Phillipps 1996). 

Blanchard et al.(1992) discuss the problem of overprc- 
dicting the small objects, and mention the Press-Schcchter 
prescription might not be the exact answer (but give no 
specific alternate prescription), blaming thus the estimate 
of the mass function. However, the same problems appear 
in numerical simulations (White et al.1987). In fact, we 
find that although the estimate of the mass function by 
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the non-linear scaling formulation does reduce the number 
of faint galaxies, our solution is mainly in implementing 
the proper effects of cooling. Blanchard et al.(1992) on 
the other hand suggest another way out, by reheating the 
universe so as to destroy the objects with the shallower 
potential. With their solution, the smaller the object the 
more efficient the destruction mechanism, and the slope of 
the galaxy luminosity function may be expected to be flat 
down to the faintest galaxies, whereas for our solution the 
slope eventually gets steep (a ~ —1.3) for the very faint 
objects. 

Kauffmann ct al.(1993) follow White & Frenk (1991) 
but focus their attention on the individual histories of 
mass condensations. They note that a fixed density thresh- 
old implies that most massive galaxies generally formed 
recently and consequently have young stellar populations 
and blue colors, at variance with what is observed. With 
our views, in total agreement with the assumption of hier- 
archical clustering, the more massive objects have a higher 
dark matter density contrast (which does not imply that 
their density within their luminous radius is larger than 
for faint galaxies) and formed at high redshift. Since then, 
they have formed stars as almost isolated systems and 
have been able to exhaust their gas, whence exhibiting 
older stellar populations. In fact, as far as star formation 
is concerned, massive galaxies have roughly the same age 
as other galaxies, that is about the age of the universe, 
since before they reached the curve C\ gas was already 
being transformed into stars, within sub-units which later 
merged to create the large galaxies we can observe today. 
In this sense, the "age" of these bright galaxies (that is 
the time available for star formation) is even slightly larger 
than the one of small galaxies, because within our model 
these massive objects are the final results of the merging 
of deep potential wells where star formation started ear- 
lier since cooling was more efficient there than in weaker 
halos. However, the main reason in our model (with the 
fact that all galaxies have roughly the same "age" ) for the 
old stellar population of luminous galaxies (as compared 
to faint ones) is that their star formation process is more 
efficient because of their high density and virial tempera- 
ture. This trend is also apparent (but to a lesser degree) 
in the model of Cole ct al.(1994). 

The study of the galaxy counts as a function of ap- 
parent magnitude, where evolution in luminosity and 
in number are intrinsically tied (Guiderdoni & Rocca- 
Volmerange 1990) is not fully satisfactory within our sim- 
ple description. It deserves a detailed study using our 
methods but including properly color evolution, and re- 
mains to be done with the needed care. Also, the fate of 
galaxies in the dense parts of clusters, when dark halos 
start to loose their individuality, would be worth a thor- 
ough consideration. Clearly, however, the obvious outcome 
of such a model where mass and luminosity evolution are 
treated consistently is the study of the galaxy populations 



at high redshift. We have here only unveiled some of the 
applications of our approach. 

Acknowledgements. We thank J. Silk for discussions during 
various stages of this work. 
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APPENDIX 

A. Characteristics of galactic halos 

Any dark matter halo of radius R can be characterized 
by its overall density contrast A or equivalently by the 
parameter x introduced in (8). We assume that these 
potential wells have a mean density profile of the form 



p{r) oc r~ 7 where 7 is the slope of the two-points correla- 
tion function £ (we shall come back to this point below). 
Since 7 ~ 1.8 is close to 2, the halo is close to isother- 
mal and the velocity dispersion a 2 oc M(R)/R oc _R 2 ~ 7 
is nearly constant, which agrees with observations. With 
such a density profile, the velocity dispersion of the dark 
matter verifies: 



_d_ 

dR 



(P« 2 ) = - 



QMp 
R 2 



which leads to 



a 2 (R) = 



27-2 



QM(R) 
R 



(Al) 



The energy equipartition kT = pm p cr 2 , where pm p is 
the mean molecular weight of the gas and m p the proton 
mass, gives for the gas temperature: 



k T{R) = 



1 



2 7 -2 



Gpm p M(R) 
R 



(A2) 



which also leads to the hydrostatic equilibrium for the gas, 
using P = ribkT = pb/(pm p ) kT, where rib and pb are the 
baryon number density and mass density, if pb oc p. Hence 
we shall assume that the gas initially follows the spatial 
distribution of the dark matter in the virialized halo, be- 
fore it starts cooling, which is consistent with numerical 
simulations (Evrard 1990). More precisely we assume: 



Pb 



where 0/, is the present ratio of the baryon density to the 
critical density. We shall take fif, = 0.04 in the case = 1, 
and fib = 0.03 in the case flo = 0.3. Both values are consis- 
tent with the bounds given by primordial nucleosynthesis 
(Walker et al. 1991). Finally, we note: 



Pb 
n b m p 



0.59 and p e = 



Pb 



1.14 



n e m p 



which correspond to halos of primordial composition (with 
an helium mass fraction Y = 0.24), where n e is the elec- 
tron number density. We note V the circular velocity of 
the galaxy at its external radius R : 



2 = QM_ = __L — 2 

R 27 - 2 P p 



(A3) 



We also define a luminous radius R c , and the circular ve- 
locity V c at R Cl by: 



R r 



1-7/2 



We shall use V c to compare our model to the observed 
Tully-Fisher relation, which relates the luminosity to the 
circular velocity at a radius ~ R c . Since 7 ~ 2 we have 
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V ~ V c (the rotation velocity is roughly constant through- 
out the halo). 

The time t m introduced in Sect. 3.1 corresponding to 
the turn-around epoch of a given halo (or to the time when 
it became non-linear) is given by the spherical collapse 
model (Peebles 1980): 



' R = A(l - cos0) 



with A 3 = QMB 2 



t = B( 



sin ( 



The time of turn-around is simply t m — irB, and the ra- 
dius of maximum expansion is R m — 2A. When the over- 
density virializes, at p t m , its radius is R v — R m /2 = A, 
and the averaged density within R v is p = M/(Aw/3 R 3 ) = 
?>ir/(Agt 2 m ). Thus we obtain: 



t r . 



3tt 



which gives the time of turn-around t m of the halo as a 
function of its density. 

The slope 7 of the dark matter halos mainly enters 
our calculations as a numerical factor of order unity in 
(A2). Hence we do not need a precise description of the 
detailed density profile of the objects we deal with. More- 
over, we explicitely consider that very large halos defined 
by A ~ 177 which correspond to clusters of galaxies (at 
low z) contain several higher-density sub-units which are 
identified as distinct galactic halos. In fact, we recognize 
directly these smaller objects without dealing (in this ar- 
ticle) with their larger host halo. 

B. Star formation model for an isolated halo 

With the prescription introduced in Sect. 3.1 we are able 
to obtain the mass function of galaxies, or the tempera- 
ture, radius functions, but we do not have the luminosity 
function yet. Indeed, the luminosity of galaxies depends 
on the processes which govern star formation, which wc 
need to take into account explicitely. Wc shall adopt a 
very simple model for the formation of stars within the 
halos we previously considered. 

B.l. Model 

For each galaxy we divide stars in two populations: a first 
class of massive stars, with a life-time T s h shorter than the 
age t of the galaxy, and a second class of small stars with 
a life-time 77 larger than t. Thus, the gas used to form 
these small stars has not been recycled in the ISM, since 
these stars are still on the main sequence, while the large 
stars consist of many successive generations which have 
continuously recycled part of their mass through super- 
novae explosions, planetary nebulae and winds. We con- 
sider explicitely the mass recycled by SNII, and by stars 
in the AGB phase, but we neglect the mass ejected by 



SNI (believed to be associated with white dwarf coales- 
cence and explosion) which is very small as compared to 
the one involved in the former processes. To include SNI 
would be straightforward, but it would not change our re- 
sults as long as we do not consider the history of Si or Fe. 
We note M g , M s h and M; D the total mass in the form of 
gas, short-lived stars and long-lived stars. Moreover, we 
consider two gaseous phases: some low-density diffuse gas 
(A5) M g h spread over the halo, and some dense gas M gc within 
the core of the galaxy in the form of clouds which turns 
into stars with a time-scale r c . Initially, we have M gc = 
and Mgh — Mi, with M^ = f2b/f2o M the mass of baryons 
in the halo (hence we assume that at the time of virializa- 
tion the baryon fraction within any halo is representative 
of the universe: there has been no prior segregation be- 
tween baryonic and dark matter). The diffuse phase M g h 
is continuously replenished by stellar winds and super- 
(A6) novae, which eject and heat part of M gc . Meanwhile, the 
diffuse gas settles in the central parts of the galaxy and 
forms dense clouds over a dynamical time-scale Td- Indeed, 
the gas ejected by supernovae or ionized by stellar winds, 
or initially hot after virialization, cools over a time-scale 
t coo i and then falls back into the center of the potential 
well over a time Td ~ t m . By definition of our halos, the 
constraints (15) ensure that t coo i < Td hence Td is the only 
one relevant time-scale. In fact we even have t coo i <C Td 
at high rcdshifts or small temperatures when p v p\, in 
this case most of the gas M g h is cold. This would be differ- 
ent for galaxies in clusters, where the gas may be spread 
all over the cluster, cooling is less efficient, and where the 
potential wells of the galactic halos are expected to be 
modified by the additional cluster gravitational energy. 
The properties of galaxy clusters, however, will be treated 
in a forthcoming paper. Finally, we can write: 



( dM gh 
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(Bl) 



where 77 <C 1 is the mass fraction of the IMF corresponding 
to short-lived stars, for 1M Q of gas which is processed into 
stars, and 77' ~ 0.1 is the fraction of mass which is locked 
within white dwarfs or neutron stars after the death of 
these massive stars, and is not recycled in the galaxy to 
form other generations of stars. We have 77 <§C 1 because 
in the case of usual stellar initial mass functions (IMF) 
most of the mass is within low mass long-lived stars. In 
a fashion similar to Kauffmann et al.(1993) we write the 
mass of gas heated and ejected by supernovae as: 



dM g 
dt 



SN 
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where e ~ 0.1 is the fraction of the energy E$n delivered 
by supernovae transmitted to the gas, and we defined: 



To 



2 e E SN fim p n SN 
3 k msN 



10 6 K 



using t)sn /m S N — 0.005 M Q X and E SN = 10 51 erg. This 
value of rjsN /tnsN is consistent with the IMF we shall use, 
and it is constrained by the observed supernovae rates. 
The effect of stellar winds can also be incorporated into 
this model through T . 

Using the fact that by definition of short-lived stars 
Tgh <fwe can make the approximation that M s h evolves 
in a quasi-static way: d *£* h ~ 0. Hence we obtain: 



M sh = V ^ M gc 



(B4) 



Then we solve numerically (Bl) to get the evolution of 
M g h, M gc and all other quantities. 

B.2. Analytic approximations 

For faint galaxies characterized by a small virial temper- 
ature T <C To the system (Bl) leads to a simple approx- 
imation, which also provides uscfull insight into the be- 
haviour of more massive galaxies, as can be checked nu- 
merically. More precisely, this approximation is valid pro- 
vided t c <C (1 + To/T) 2 Td- In this case, a quasi-stationary 
regime sets in very quickly, where the galactic evolution 
is regulated by supernovae (and stellar winds) and the 
dynamical time-scale: the sink term for the star-forming 
gas Mg C , corresponding to matter ejected by supernovae 
(which involves the star formation rate and the super- 
novae efficiency parameterized through To), balances the 
source term due to the infall of gas from the extended halo 
Mgh (which involves the dynamical time-scale). Hence the 
global star formation rate is governed by the interplay of 
the supernovae efficiency and the dynamical time-scale. 
Then, M gc evolves in a quasi-static way and follows closely 
the mass of the reservoir M g h- 



Mgc = [1 + T /T - r/(l - r/)]" 1 r c /T d M, 



gh 



(B5) 



Note that the condition of validity of this approximation 
implies that M gc <§C M g h- Then, since M g ~ M g h we ob- 
tain 



M a = M g0 e-W T ° 



where t g is the age of the galaxy and we defined: 



1 + Tp/T - 7/(1 - rf) 

TO = : 7Z J\ T d 

l-?7(l-r7') 



1 • y I r,, 



(B6) 



(B7) 



Thus, as we explained above, the galactic evolution is gov- 
erned by T /T and Td, the time-scale r c does not appear in 



the global gas mass, or stellar content. The instantaneous 
star formation rate is 



dM s _ M g 
(B3) ~dT- — 



(B8) 



since M s + M g = M g0 , hence dM s /dt = -dM g /dt (we 
neglect the mass in the form of short-lived stars). The 
mass within massive short-lived stars is: 



M s h = , 1 jt — M g 

1 - 7/(1 - 77') T 



(B9) 



while the mass in the form of small long-lived stars is: 
I-77 



M lo 



1-77(1-77') 



- (M s0 - M g ) 



(B10) 



The mass locked in star remnants, white dwarfs or neutron 
stars, is 



M r = -f^- M lo 

1-77 



(Bll) 



C. Redshift evolution and merging of galaxies 

C.l. Model 

In the previous section we considered non-evolving halos, 
in the sense that the total mass of the galaxy, and its time- 
scales t c and Td, remained constant with time. However, as 
we explained in Sect. 3. 2, the characteristics of the galactic 
halos we consider evolve with time, following the curve C 
(see Fig. 2). Thus, we now have to model the evolution 
of the matter located on C v through its merging history. 
Similarly to White & Frenk (1991), we can write for the 
comoving stellar content of the universe at time t: 

r* r°° dM 

7]{x,t)M s (x,t) = dt' dx' —~^{x' ,t') 
Jo Jo dt 

x n(x',t') P{x,t\x,t') (CI) 

where P(x,t\x' ,t') dx is the probability (which we do not 
know) that matter which was in a halo of parameter x' at 
time t' will be part of a halo of parameter x - x + dx at 
t, and dM s /dt(x',t') is the star formation rate of the cor- 
responding halo. We can define halos by their parameter 
x for both PS and non-linear scaling approaches because 
we noticed in Sect. 2 that the usual linear parameter v is a 
function of x only, see (13). If individual halos are stable 
once in the non-linear regime (when A > A c and £ ^> 1), 
their parameter x remains constant with time, but they 
accrete some mass as their virialization radius gets larger 
with time and they may join together to form a larger 
halo. Hence, we shall make the approximation that the 
mass which is within halos x — x + dx at t was within 
halos of the same range in x at all earlier times. This also 
satisfies the constraint that mass be conserved since the 
mass fraction within halos x — x + dx does not depend on 
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time, and is a function of the sole variable x, see (10). In 
fact, the conservation of mass implies to use a distribu- 
tion in x, or to take x constant which should be a good 
approximation if this distribution is peaked around a par- 
ticular value. This also ensures that we do not mix the 
mass of non-evolving galaxies located on Ca with the ha- 
los on C v , and that the mass at a given time t comes from 
older less massive halos characterized by a slightly weaker 
potential well since for a fixed x the virial temperature T 
is smaller at higher redshifts. This is quite natural since as 
time goes on potential wells merge to form deeper ones. In 
other words, we neglect the scatter in all possible "merger 
trees" and simply use a "mean history" defined by com- 
patibility requirements. This analysis also applies to the 
PS approach, which would only give a different function 
h(x), hence a different mass fraction. 

In fact, when t' -C t one expects that the correlation 
x' — x disappears, or at least gets weaker, but this is not 
a problem if the stellar content of a given object at t is 
dominated by its recent star formation history, which is 
the case. Moreover, within the framework of our model for 
star formation small galaxies are regulated by the inter- 
play between supernovae (which eject gas) and the infall 
of gas from the outer parts of the halo and their lumi- 
nosity is dominated by recently formed stars which are 
more numerous. As a consequence the details of their pre- 
vious stellar history are not very important. On the other 
hand, very massive and bright galaxies are not affected 
by the supernovae feedback mechanism (since their po- 
tential well is sufficiently deep to retain the gas very effi- 
ciently) so that their past history matters. Besides, at low 
redshifts their star formation rate begins to decline sig- 
nificantly as they exhaust their gas content. This means 
that their luminosity and stellar properties are governed 
by old stars which formed during ealier and more active 
periods. However, in our astrophysical model we assume 
that these galaxies (located on the cooling curve Ca) do 
not evolve significantly any more so that our approxima- 
tion x' = x becomes correct for these objects. Thus, our 
approximation is consistent with our model and it should 
provide a reasonably good description. Moreover, it allows 
one to get simple analytic insights into the global galaxy 
formation process, which clearly show the general trends 
implied by any such model based on hierarchical structure 
formation supplemented by cooling constraints. Note that 
using a different galaxy formation model Kauffmann et 
al.(1998) found that their results were not very sensitive 
on the details of the merging trees of their halos (they 
obtained similar results with N-body simulations and an 
extended Press- Schechter theory for the properties of indi- 
vidual galaxies although the merger trees differ in detail). 
Using (9) we obtain 



dt ' { M 



1 dM, 



dt 



(C2) 



In fact, this approach simply means that we can still use 
the sytem (Bl) to get the proportion of the mass of gas 
which is converted into stars in a galaxy, but the time- 
scales r c and Td, and the virial temperature T, are now 
functions of time. 

C.2. Analytical solutions 

Although in practice we compute numerically the solu- 
tion of the system (Bl), we present now the case of the 
simplified quasi-stationary regime, corresponding to small 
galaxies, to give a clear illustration of the effects of this 
additional time-dependence. Moreover, since the temper- 
ature T decreases at higher redshifts, for a fixed x, all 
galaxies follow this regime when they are young. 

Since M s + M g = M b = fi b /O M, and dM s /dt = 
M 9 /to in this case, we obtain: 
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M, 
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(C3) 



Thus we have an equation similar to the one describing 
halos located on Ca, but now the global star formation 
time-scale To depends on t' . The density p v of the halos 
located on C„, which virialize at the time t V i r — p t m , is 
according to the spherical model: 



Pv = 



p 2 3ir 



(C4) 



Hence p v oc t~ 2 , whatever the value of flo. If clustering 
is stable, for £ ^> 200, the quantity p£(R,t) is constant, 
hence £ oc i?~ 7 p ~ . Thus, for a constant parameter x = 
Pv/(pO w e obtain for halos on C v : 



R oc p- 1 ^ and T oc p 1 ^ 2 ^ 



(C5) 



As a consequence, since To oc p~ 

-i/2 T -i for halog guch that 

T < T , see (B7), we get: 



r (i) oc p- z ' 2+2 ^ oc t 3 - 4 ^ along C v at fixed , 



(C6) 



Hence, if we note tn the age of the universe at the redshift 
we consider, and To the star formation time-scale at this 
date, we can write: 



M, 



so 



H dt { M, 



r H dt 
Jo T o 



M, 



</0 



2(1--) 
3(3 + n) 



(C7) 



using 7 = 3(3 + n)/(5 + n). We can note that this integral 
converges for t —-> 0, and T (which measures the depth of 
the potential well) decreases at higher redshifts, provided 
— 3 < n < 1, which corresponds to the range of interest 
where hierarchical clustering is valid. Hence our analysis 
applies to all relevant power-spectra P(k). Thus we obtain 
for the gas mass fraction at any time t: 



M a 



M, 



= exp 



90 



3(3 + n) t_H_ (t_ 
2(1 -n) t \t H 



2(l-n)/(9+3n) 
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In the case n = — 2 we have for the gas mass fraction at 
time iff, when we calculate r , 

Ma_ = e -W(2r ) on Cv with t<Tq (CQ) 

MgO 

Thus, we obtain a relation similar to (B6), with an "ef- 
fective age" t g for the galaxy given by t g = t H /2. Hence 
the time-evolution of T (t) reduced this age by a factor 
1/2, as compared to (B6), since for n = — 2 the time- 
scale T (t) increases at high redshifts, which leads to a 
less efficient star formation, because the temperature T 
decreases. This would not be the case for large n, where 
the dominant effect would be the increase of the density. 
Note however that the variation with n of the factor 1/2 
can simply be incorporated into the parameters (3 C and 
(id- Galaxies located on C v with a high virial temperature 
T > T have a slightly more complicated history. There is 
a first stage, at small times, where T < T and they follow 
the behaviour we have just described, and a second stage 
where T > T and r (t) oc t. Hence we divide the integral 
of (C3) in two parts, and we get a similar result, with an 
effective age: t g = [1 + ln(T/T )] t H /2. In fact, the rela- 
tion (B6) does not apply to these massive galaxies, but it 
still gives a good estimate of their evolution, and this sim- 
ple result shows that the effective age is increased because 
the temperature remains constant for some time, as could 
be expected, but only by a logarithmic term. For massive 
galaxies located on Ca we have To = constant down to 
the time t vir when A = A c and we switch onto C v . At 
small redshifts z < z V i r the fact that To is constant means 
that the gas mass fraction follows a simple exponential 
decline M g cx exp(— t/r ), usually different from the evo- 
lution M g cx exp[-0/To) 2(1 ~™ )/(9+3n) ] along C v . The time 
of virialization is simply t v i r = p t m and the star forma- 
tion time-scale at this date is to- Hence these galaxies have 
an effective age t g which is the sum of their effective age 
at t v i r , which we obtained above, and of the time which 
has elapsed since this date: tn — P t m . Thus, to sum up, 
for all galaxies the quasi-stationary approximation gives 
again the relations of the previous section, (B6) to (Bll), 
with an effective age given by: 

( T < T : t g = 1/2 p t m + (t H -p t m ) 

[t > To : t g = 1/2 [1 + In(T/T )] p t m + (t H -pt m ) 

Note that halos on C v satisfy tn = P t m by definition. 
Thus, the gas mass fraction within halos at time tn is 
given by a relation of the form M g /M g o — cxp(— tn/ro), 
whatever the precise dependence on z of the star formation 
time-scale, provided the integral of the right-hand side of 
(C3) converges (star formation does not occur as a sudden 
burst at high redshifts). Moreover, the variation with n of 
the factor 1/2 in the age of the galaxy, when it is on C v , 
is simply incorporated into the parameter (3d which enters 
the definition of t , see (B7). Thus, as far as star forma- 
tion processes are concerned, all galaxies have roughly the 
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same age (of the order of the age of the universe), even 
if they did not exist as distinct entities over this whole 
period of time: they continuously accreted some mass or 
merged with neighbours. 

We can see that the time evolution of these virialized 
halos does not depend on flo- This is due to the fact that 
the internal dynamics of overdensities given by the spher- 
ical model does not depend on the background universe. 
However, the number of such halos will naturally vary 
with the cosmological parameters. Besides, the rcdshift 
evolution depends on flo, in parallel to the relation time- 
redshift. For instance, if O = 1 we have (1 + z) oc t~ 2 / 3 
while (1 + z) oc i" 1 if SI ~ and A = 0. Thus, in the 
case n = —2 we have for halos on C v defined by a fixed 
parameter x: 

'{1 = 1 ■ Ml. = e -V(2ro) (l+*r 3 

MgO 

(CIO) 

Q « 1 , A = o : ¥°- = e -V(2n,) (i+z)- 2 

l Mg 

where t = tn{z = 0) is the present age of the universe. In 
the case of a low-density universe the redshift evolution is 
slower, because of the faster expansion which implies that 
the same multiplying factor in redshift corresponds to a 
smaller factor in time. 

We can note that for halos on C v , the time-scale which 
governs their merging history is r^. Hence, if star forma- 
tion is also enhanced by gravitational interactions with 
surroundings and merging with other halos, the natural 
time-scale is again which gives another justification for 
our star formation time-scale t c cx tj. The gas of halos 
which undergo this succession of mergings can be reheated 
to the virial temperature by the energy released during 
these violent encounters. However, for halos located on C v 
we have by definition t coo i -C r<j, hence we can neglect 
these successive phases of reheating. 

The evolution equation (C3) we obtained is interesting 
as it readily shows the influence of the parameterization 
adopted for the star formation time-scale r on the history 
of galaxies. Indeed, since p v cx t ^, we can see that the in- 
tegral on the right-hand side diverges for t' — ► if To is a 
strong power-law of the density: To cx p~ a with a > 1/2. 
In fact, there is a cutoff at high redshifts because for a fixed 
parameter x the temperature gets smaller in the past, and 
when T < 10 4 K cooling is very inefficient so that star 
formation is suppressed. However, this large decrease of T 
by a factor 100 (from <~ 10 6 K to 10 4 K) means that the 
redshift of cutoff z cut corresponding to present galaxies is 
rather high: z cut ~ 4. Thus, such a model for r would 
imply that most stars formed at z ~ 4 and that the global 
star formation rate has been very low ever since. This is 
certainly inconsistent with observations, which show that 



32 



Valageas & Schaeffer: The Mass and Luminosity Functions of Galaxies and their Evolution 



the star formation rate has not experienced such a dra- 
matic decline since z ~ 4 and that star formation is still 
active in the present universe. If To varies weakly with p: 
a < 1/2, the integral converges and star formation is dom- 
inated by the latest epochs, and increasingly so for smaller 
a, until all the gas is converted into stars (at some time 
in the future). The case a = 1/2, which corresponds to 
our prescription, is intermediary as the divergence would 
only be logarithmic, which is not a problem because of 
the cutoff. In fact, the additional temperature dependence 
(1 + Xb/T) introduces another rcdshift factor which makes 
the integral converge. This strongly suggests that the star 
formation time-scale To should vary as t oc p^ 1 ^ 2 , at least 
at high redshift, (with a possible additional dependence on 
T), independently of the arguments presented in the pre- 
vious section, so that star formation is neither dominated 
by a sudden burst at high redshift when most of the mass 
reaches T > 10 4 K, nor by the very recent epochs, since 
these both cases are inconsistent with observations. Note 
that it is correct to consider the quasi-stationary regime in 
this analysis since we are interested in the earliest stages 
of galactic evolution when T <C T . 

C.3. Total mass in stars and global star formation rate 

We can also consider the evolution with rcdshift of the 
total mass in stars per comoving Mpc 3 . The mass of 
baryons within galaxies is J Qb/Qo M n(M) dM/M. Us- 
ing M b = M g0 = M g + M s we obtain: 



since To oc p~ 1,/2 , see (B7), and p oc (1 + z) 3 for most galax- 
ies which are close to C v , we obtain dp s /dt oc (1 + z) 3 / 2 as 
long as T > Tq. However, at high rcdshifts the comoving 
star formation rate gets smaller because the mass con- 
tained in deep potential wells starts to decrease (influence 
of the cutoff at T = 10 4 K), and more importantly because 
as the virial temperature declines the star formation time- 
scale To starts to increase (in the case n = —2), see (B7). 
Naturally, if we neglect the very small apparent mass loss 
at T < 10 4 K, we have: 



dp s d 

~* = dt mz) 



Pc(z)] 



(C14) 



by construction, as implied by (C2) (and checked numer- 
ically). 

We can also look at the evolution with rcdshift of 
the star formation rate dM s /dt of individual halos de- 
fined by a fixed temperature T. Note that two such ha- 
los defined by the same temperature at different rcdshifts 
are not necessarily formed by the same matter. Along 
C v , we have p oc (1 + z) 3 , hence t oc (1 + z)~ 3 / 2 , and 
M oc (1 + z) -3 / 2 . For ft ~ 1 the age of the universe scales 
as tn oc (1 + z)~ 3 / 2 , and exp(— t H /2r ) ~ 1, hence we 
obtain dM s /dt = 1/t M g0 exp(— t/2r ) ~ constant with 
time. In a similar fashion, for halos located on C\ we also 
get dMg/dt ~ constant as long as t < t . Thus, the star 
formation rate is roughly constant with time for these ob- 
jects, which is consistent with observations. 



n 



Ms 



M, 



p(M) 



dM 
~M 



(Cll) 



where Cl s (z) = p s (z) / p c (z) is the stellar density parameter 
and p(M)dM / M is the mass fraction in galaxies of mass 
between M and M + dM, given by (3) and (10) for the PS 
and non-linear scaling prescriptions. Naturally, we always 
have n s (z) < f2&/Oo £l(z),asis clearly seen in (Cll), since 
J p(M)dM/M — 1. In fact, there is also a low temperature 
cutoff at T = 10 4 K, but up to z ~ 4 this only involves a 
very small fraction of the total mass. 

We obtain the comoving star formation rate in the 
same way. We have dM s /dt = [1 — r/(l — ?]')] M 9C /t c , for 
long-lived stars, hence: 



dp s 
dt 



77(1 



/m m qc ,*,s dM 



(C12) 



For the quasi-stationary regime we get for the derivative 
of the stellar comoving mass density: 



dp s _ Mb f J_ 
dt O J t 



f /-r . . dM 
- t/T0 KM) — 



(C13) 



It first increases with redshift because the gas content of 
bright galaxies gets higher (depletion term exp(— t/ro)) 
and the star formation time-scale tq decreases. Indeed, 



D. Stellar properties of galactic halos 

We can note that the mass in the form of short-lived stars 
is always much smaller than the mass contained in long- 
lived stars, and increasingly so for well evolved galaxies 
(when t » t ) as could be expected. Indeed these galaxies 
have already consumed most of their gas content, hence 
their present star formation rate is relatively small (as 
compared to their past) and their stellar population is 
dominated by all stars which were formed all along the 
galaxy history. Thus, for the quasi-stationary approxima- 
tion we obtain: 



Ms 



Mio 
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Tsh 
TQ 



exp 



< 1 



(Dl) 



where we used i)« 1, since in the case of usual stellar 
initial mass functions (IMF) most of the mass is within 
small long-lived stars, and we have r s h *C t. Thus, lumi- 
nous galaxies, which correspond to large circular veloci- 
ties (as is observed through the Tully-Fisher relation) and 
high densities (because of the cooling constraint, see the 
curve Ca in Fig. 2), hence to a small global star formation 
time-scale t , have consumed most of their gas and will 
be redder than faint galaxies. This is an important success 
of our model as this trend is actually observed (Lilly et 
al.1991, Metcalfe et al.1991), but usually difficult to get 
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by common models. We can also notice that the mass in 
the form of white dwarfs or neutron stars is always very 
small as compared to the total stellar mass from (Bll). 
Indeed we get (with 77 ~ 0.1 and 77' ~ 0.2): 



M r 
M7 



1 



0.02 



V 



(D2) 



as is the case in any reasonable stellar evolution model. We 
can check that mass is conserved with time in the sense 
that: 



M g + M io + M r = M g0 



(D3) 



In fact we have a slight excess of mass M s h, which is 
negligible since we noticed above that M s h <§C M; G for 
instance. This is due to our quasi-static approximation 
^f- ~ 0. Finally, we also consider that for 1M Q of 
gas converted into stars a fraction rjd goes into "invisi- 
ble" compact objects, such as brown dwarfs, which we 
included among the long-lived stars. Since observations 
seem to show that this mass fraction is rather small we 
choose rjd = 0.1 in the numerical calculations (in fact we 
could as well use % = or 77^ = 0.2, since this parame- 
ter has almost no influence as long as it remains small). 
Then, the mass M s of luminous stars is in our model 
M s = M sh + (1 - 77 - r? d )/(l - ??) M to ~ M to , since we 
have to remove the part of M; c formed by dark objects, 
but this fraction is negligible. Besides, 77 <C 1, 77^ <C 1 
and we have already noticed that M s h <C Mi a . Thus, we 
obtain: 
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(D4) 



At small times when the gas content of the galaxy is still 
important we have: 



Mg_ 



TO 
t 



for t < t 



(D5) 



It is clear on this expression, which is valid for any model 
of star formation and does not rely on the quasi-stationary 
approximation, that the observed ratio M g /M s of the 
Milky Way and its age give directly its star formation 
time-scale To (whence the parameter (3d)- 

Thus we now have attached a peculiar stellar content 
to each halo, or galaxy. To get the luminosity of such a 
galaxy, we only need to precise the luminosity of its stars. 
We note L s h and L/ Q the global luminosity of short-lived 
and long-lived stars per unit mass, so that the luminosity 
L of the galaxy is: 



L = L sh M sh + L lo M lo 



(D6) 



Now we have to precise the values of the quantities 
77, 77', T s h, L s h and Li a . We shall derive them from the ini- 
tial mass function (IMF) of stars and mass luminosity and 



mass life-time relations. For 1M© of matter converted into 
stars, the number of stars dM* formed in the mass range 
to — m + dm is 



d/V* = <t>{m) dm with <p(m) — am 



(D7) 



where m is the star mass in units of 1M Q and a the nor- 
malization constant. We use x = 1.35 for m < 1 and 
x = 1.6 for m > 1, which is similar to the IMFs given 
by Salpctcr (1955) and Scalo (1986). This applies to stel- 
lar masses between to_ = 0.1M Q and m + = 125M Q . 
We could change somewhat this IMF (for instance choose 
x = 1.8 for m > 1) without significant variation in our 
results. Moreover, a fraction r\d may go into "invisible" 
compact objects, such as brown dwarfs, so we have 



i = m + 
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J m 



m <t>(m) dm 
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which defines the normalisation of 4>{m). We note to* the 
mass which separates the two classes of stars we intro- 
duced above. Then, we obtain: 



77 = / to <j){m) dm 



and 



(1 - 77) = r]d+ / m<t>{ 
J Tu- 



rn) dm 



(D9) 



(D10) 



The mass M; G is formed of luminous stars, with m_ < 
to < to*, and of dark objects. Next we use the mass lumi- 
nosity and the mass life-time relations: 

L = Lq m v and n = 10 10 ^ = 10 10 m 1 ^ years 

Mq Li 

with v = 3.3, which is consistent with the mean observed 
mass — B band luminosity relation for stars on the main 
sequence. Hence we have: 



i] T s h = m Ti(m) (j){m) dm 



(Dll) 



and 



V T s h L s h = / 17 (to) L(m) 4>(m) dm (D12) 

J 

Similarly, we get: 

(1 - 7]) Lio = / L(m) 0(m) dm (D13) 

Using (D7) and the previous mass-luminosity relation we 
can see that the luminosity of the galaxy is dominated by 
the contribution of stars of mass close to to*. Finally we 
use 77' = 0.2, and the mass m* which separates the two 
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classes of stars we introduced above is chosen to be given 
by: 



r* = T;(m*) = 0.5 t H 



so that T s h is smaller than the age of the galaxies, which is 
close to the age of the universe at the time of interest. For 
instance, in the case Q = 1 a galaxy like the Milky Way, 
that is with a circular velocity V c = 220 km/s, corresponds 
in the present universe to to = tji{z = 0) = 11 10 9 years, 

= 1.3 , Tgh/to = 0.11 and r\ = 0.25. 

Thus we can attach a luminosity to each galaxy, which 
enables us to get the luminosity function of galaxies from 
the mass function: 



r,{L) - = v(M) — 



(D15) 



We can also use this model to obtain the supernovae 
rate in galaxies. Thus, we assume that stars more massive 
than m,i = 6 M Q will explode as supernovae after they 
leave the main sequence. Naturally, they must also be part 
of our class of short-lived stars to explode during the life- 
time of the galaxy. The mass in the form of these massive 
stars is: 



M S n = Vsn M, 
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(D16) 



Thus, in the case of the quasi-stationary regime we obtain 
for the number of such stars at any time: 
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and the supernova rate is: 
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(D18) 



(D19) 



In the case fl = 1, we obtain for a galaxy similar to the 
Milky Way a supernovae rate of 2.5 explosions per century, 
which is close to the observed value for type II supernovae. 

E. Metallicity 

Finally, we can derive the metallicity from our model (Bl). 
The metallicity here is understood as being the abundance 
of Oxygen, or any other element that is not significantly 
produced in SNI which are not included in our model. 
We define Zh and Z c as the fractions of metals within 
the diffuse gas M g h and within the gas in the core of the 



galaxy M gc . We note Z s the mean stellar metallicity. Thus, 
we obtain: 



( dZ h 
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dZ s 
dt 



[1 - 7,(1 - 7/)] 



Mgc 
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where y is the yield. The fractions of metals in both gas 
phases M gh and M gc vary because of the exchange of mat- 
ter between these two components. In addition, the central 
gas M gc is enriched by stellar ejecta. The mass of met- 
als in stars increases as metals are incorporated from the 
star- forming gas M gc . In the case of the quasi-stationary 
regime, we obtain for the gas in the diffuse phase: 



Z h = 
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to 



1 + T/T T 



while we get for the dense central phase: 
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where t g is the age of the galaxy. Since the halo is not en- 
riched directly, but through the mass loss of the galactic 
core, its metallicity remains for a long time much smaller 
than the one attached to this central component which 
receives the stellar ejecta. The second term of Z c , which is 
constant in time, corresponds to the fact that very quickly 
a "quasi equilibrium" regime sets in where the gain of 
metals within this central phase M gc from stellar ejecta 
balances the loss due to the exchange of matter with the 
diffuse phase, which replaces some gas with the metallicity 
Z c by some gas falling from the halo with the metallicity 
Zh <C Z c . Since this exchange of gas between both com- 
ponents Mg C and M g h is driven by supernovae, or stellar 
winds, whose importance is parameterized by Tq/T, see 
(B2), the equilibrium metallicity is smaller for small tem- 
peratures (weak potential wells) where these flows of mat- 
ter are more important. Then Zh increases linearly with 
time as M g h receives some metals from M gc at a constant 
rate (the term in Zh corresponds to the second term in 
Z c , simply multiplied by a temperature-dependent factor 
and time). Naturally, when Zh reaches the "stationary" 
value described above for Z c , this "equilibrium" regime 
stops, both phases have the same metallicity which in- 
creases linearly with time as the enrichment process goes 
on. The stellar metallicity is: 
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1 + T/T 



exp(i g /T ) - 1 



+ 
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1 + T /T 



(E4) 



The second term in Z s , which is constant in time, cor- 
responds to the "stationary" regime we described for Z c - 
Since stars form from the dense phase M gc , their metal- 
licity quickly reaches the surrounding gas metallicity Z c 
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which is constant with time. Later, when most baryons 
are incorporated into stars (t g > To), the stellar metallic- 
ity increases to reach y if the previous equilibrium value 
was smaller than y. Indeed, at very long times when all the 
gas is converted into stars the mean stellar metallicity Z s 
must be equal to the yield by definition, since our system 
is closed (there is no global mass loss nor gain although 
the mass of individual galaxies evolves with time). 

F. Role of parameters and scalings 

The model we constructed in the previous sections has 
only two main specific parameters: (sq) and (3d, which 
determine respectively the size (and mass) of galaxies 
and their star formation rate. For massive galaxies, with 
T > Tq, there is also a dependence on (3 C , while for faint 
galaxies there is a dependence on To, but this is not crit- 
ical. Naturally, it also depends on the usual cosmological 
parameters tto, Qb, h, and on the power-spectrum P(k) 
for the multiplicity functions. Now we shall see how we can 
get the value of (sq) and (3d, for a given set of cosmological 
parameters and a fixed To. 

Let us consider a galaxy like the Milky Way located on 
the curve C\ (see Fig. 2), defined by a fixed temperature T 
or circular velocity V. From the relation t coo i = q t m we 
obtain p oc (sg)~ 2 ilg fl b 2 , and then all the characteristics 
of this galaxy. Thus, we get: 



^ ex f3 d (sq)n- 1/2 n h A^h 



and 



L io M io cx (sq) 



(Fl) 



(F2) 



Besides, since M s h <C Mi Q we have L <~ T; D Mi Q , and in- 
creasingly so for massive galaxies like the Milky Way as 
we noticed earlier (in fact, for faint galaxies which have 
not exhausted their gas content L s h M s h ~ Li a Mi a while 
for bright galaxies which have already consumed most of 
their gas content L s h M s h <C Li a M/ Q ). Hence, a larger 
(sq) leads to a higher gas/star mass ratio and a larger lu- 
minosity, which is quite natural since it means a less con- 
straining cooling constraint, whence a broader, more mas- 
sive and lower density halo. Similarly, a larger (3d leads 
to a higher gas/star mass ratio, since it means a longer 
star formation time-scale, whence fewer stars. It does not 
influence the luminosity because the mass in the form of 
long-lived stars is of the order of the initial mass of gas for 
these dense galaxies. As a consequence, the observed lu- 
minosity and gas/star mass ratio of the Milky Way, which 
corresponds to a circular velocity of 220 km/s and a tem- 
perature T = 3.6 10 6 K, give the value of (sq) and (3d, as a 
function of the cosmological parameters. Thus we obtain: 

{sq) oc and (3 d oc % ^A" 1 / 2 ^ 1 (F3) 



Now, we can use these values of (sq) and (3d to get the 
variations of all physical characteristics of the galaxies 
we considered in our model with the cosmological pa- 
rameters. Some of these scalings are shown in Tab.l. 
However, these relations are quite general. For instance, 
the fact that M gc /M s ~ 0.4 for the Milky Way gives 
M s ~ 0.7 0(,/f2o M if the initial ratio of baryons in 
the halo is representative of the Universe, there has been 
no loss of baryons since the formation of the galaxy 
and most of the gas is in the disk. Then, the luminos- 
ity Lq of the Galaxy, and its ratio (M s /L)g, lead to 
M ~ Oo/^{, 1.3 (M s /L)g Lq oc Clo/flf,. In this way we 
obtain the mass and the radius of the halo from which the 
baryons which constitute the Milky Way came, whence 
the position of the curve C\ and the value of the product 
(sq). For fl = 1, Q b = 0.04, we get M = 2.6 10 12 M . In 
Tab.l we consider two types of galaxies: i) massive galax- 
ies like the Milky Way located on the cooling curve C\, 
which have a high luminosity and have already consumed 
most of their initial gas, and ii) small galaxies located on 
C v (hence characterized by the density contrast A c ), which 
arc faint and have a small stellar content. 



Table 1. Scalings 
Massive gj 



Small galaxies 



p oc (Q. b /Q. ) 2 
R oc Q, /Q, b 

m oc n /fit 

p /M oc 0, b h 2 



p oc (A c flo h 2 ) 
Roc (A c Q h 2 )- 1 / 2 
M cc (A c fi h 2 )- 1/2 
po/M oc Al /2 n* /2 h 3 



Thus, the values of the cosmological parameters im- 
ply well-defined galaxy characteristics. For instance, for 
a given £lo, a larger £l b (that is a higher baryon frac- 
tion) leads to a smaller radius R (since the total baryon 
mass, linked to the luminosity, must not vary), a smaller 
mass M, whence a larger number of objects (measured 
by po/M). Since the radius of a galaxy similar to the 
Milky Way should be larger than 60 kpc we get a higher 
limit on while the observed luminosity function, which 
gives the number density of galaxies, sets a lower limit. 
These two limits could have been incompatible, since for 
Oo = 1 there is no freedom in the choice of the func- 
tion h(x) which determines the multiplicity function in 
the non-linear scaling approach. Thus it appears to be 
a remarkable success that for fio = 1 it is possible to 
find a baryon fraction (in our case = 0.04) which is 
compatible with the three constraints provided by 1) the 
nucleosynthesis predictions, 2) the galaxy masses and 3) 
their luminosity function. Note however that for this latter 
constraint we still have the choice of the power-spectrum 
P(k), that is its local index n and its amplitude parame- 
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terized by as- Nevertheless, these two parameters are also 
constrained by observations and we find that the choice 
<j% = 0.5 and n ~ —2 or a CDM power-spectrum (which 
is in good agreement with observations) provides satisfac- 
tory results. In fact, we could also choose a lower value 
for fife, but this would create some problems for Lyman-a 
clouds (Valageas, Schaeffer & Silk 1998). In the case of a 
low-density universe U = 0.3, A = 0, the conditions 1) 
and 2) imply Q b ~ 0.03. 

G. Approximate power-law regimes 

From the locus of galaxy formation implied by virialization 
and cooling we may distinguish using the quasi-stationary 
approximation three regimes characterized by a specific 
power- law behaviour: 

1) Very faint galaxies, located on C v : A ~ 200 and 
T < 10 5 ' 5 K, with a constant density contrast. 

2) Faint galaxies, located on C\ : R ~ 120 kpc and 
10 5 ' 5 K < T < 10 6 ' 5 K, with a nearly constant radius. 

3) Bright galaxies located oiiCa : R — 100 kpc and 
T > 10 6 - 5 K, which nearly have exhausted all their gas. 

G. 1 . Star formation 

1) Very faint galaxies 

Using (A3), (B7) and (D5) (since these galaxies have 
only consumed a very small fraction of their gas content) 
we get: 



M cx V 3 and ^ cx V~ 2 



(Gl) 



Moreover, since these galaxies have a small stellar content 
M s < M g we have M g ~ M b cx M, hence: 

M, cx V 5 



2) Faint galaxies 
We obtain: 

M x V 2 , ^ cx V~ 3 and M s cx V 5 



(G2) 



(G3) 



3) Bright galaxies 

Since these galaxies have lost most of their gas, which 
has already formed stars, we now have M s ~ M b cx M, 
and using tq cx p -1 / 2 cx V^ 1 we get: 



Mcx V 2 



Mg 

M s 



cx exp{-V/V ) 



(G4) 



G.2. Mass/Light ratio and Tully-Fisher relation 



For a constant star-mass/luminosity ratio M s /L we obtain 
for the 3 regimes we described above the scaling relations: 
1) Very faint galaxies: 



2) Faint galaxies: 
LxV 5 and L cx M 5 / 2 

3) Bright galaxies: 
LocV 2 and L cx M 

G.5. Metallicity 



(G6) 



(G7) 



1) Very faint galaxies: 

Since we noticed earlier that L cx V 5 , we obtain: 



Z h x V 2 x L 0,4 and Z s ~ Z c x V 2 x L 

2) Faint galaxies: 

We still have L cx V 5 , so we get: 

Z h xV 3 xL ' 6 and Z s ~ Z c x V 2 cx L 



0.4 



(G8) 



(G9) 



The increase of the slope of Zh, as compared to the previ- 
ous regime, is due to the fact that the evolution time- 
scale To gets smaller for more luminous galaxies (their 
density increases, which implies that their dynamical time 
decreases), and for a constant galactic age (of the order 
of the age of the universe) it means that the halo is more 
evolved, hence more enriched. 
3) Bright galaxies: 

These galaxies have already converted most of their 
initial gas content into stars, which implies that Z s ~ y. 
Most of the gas is in the dense component M gc , since 
supernovae are not very efficient and to > tj, hence Zh — 
Z c and we recover the usual one component closed-box 
model, with Z = y ln(M ff0 /M ff ). Thus, in our case we 
obtain Z c ~y t/r . Since L x V 2 we have: 



Z r xVcxL ' 5 and Z, 



G.4- Slope of the luminosity function 



(G10) 



LozV 5 and L cx M 5 / 3 



Using the PS approximation, we get for an n = —2 initial 
spectrum (that gives results quite close to those obtained 
with CDM initial conditions) the following behaviour. 

1) Very faint galaxies: 

r){L) x L ( "- 3)/10 hence r)(L) cx L~ 5 

2) Faint galaxies: 

r){L) x L {n -V' lb hence r)(L) x L~ 0,2 

3) Bright galaxies: 

V (L) X L (n-D/e e -(L/L 3 r^ and v{L) K L -0.5 e -L/L 3 

Using the non-linear scaling approximation, we write 
h(x) x x^~ 2 for x <C x*, where uj = 0.4, and h(x) x 
jj-Wj-l e -x/x, £ or x ^ x>t w ith w s = —1.4. 



(G5) 
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1) Very faint galaxies: 

r](L) oc ^(t" -3 )/ 5 hence r](L) oc L~ 0A6 

2) Faint galaxies: 

r]{L) oc L 2 ^ -1 )/ 5 hence ry(L) oc L" ' 24 

3) Bright galaxies: 

r)(L) oc L"*e- L I L * hence ry(_L) oc L~ 1 - 4 e~ L / i * 
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